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Abstract. Stochastic point processes relevant to the theory of long-range aperiodic order 
are considered that display diffraction spectra of mixed type, with special emphasis on explic- 
itly computable cases together with a unified approach of reasonable generality. The latter 
is based on the classical theory of point processes and the Palm distribution. Several pairs 
of autocorrelation and diffraction measures are discussed which show a duality structure 
analogous to that of the Poisson summation formula for lattice Dirac combs. 



1. Introduction 

The discoveries of quasicrystals [S3], aperiodic tilings [13 [3S], and complex metallic alloys 
[60j have greatly increased our awareness that there is a substantial difference between the 
notions of periodicity and long-range order. Although pinning an exact definition to the con- 
cept of long-range order is not yet possible (nor perhaps desirable at this intermediate stage, 
compare the discussion in [56]), there is still some general agreement that the appearance of 
a substantial point-like component in the diffraction of a structure is a strong, though not a 
necessary, indicator of the phenomenon. 

Mathematically, the diffraction - say of a point set A in 1R 3 - is the measure 7 on M 3 
which is the Fourier transform of the volume averaged autocorrelation 7 of A (or, more 
precisely, of its Dirac comb 5a = J2xeA^)- Over the past 20 years or so, considerable 
effort has been put into understanding the mathematics of diffraction, especially conditions 
under which A is pure point diffractive, in the sense that 7 is a pure point measure, compare 
[aSHISlESlElElESllMlHa]. At this point in time, we have a good collection of models for 
producing pure point diffraction, particularly the cut and project sets (or model sets). Under 
certain types of discreteness conditions, one can even go as far as to say that these types of 
sets essentially characterise the pure point diffractive point sets [6]. 

But real life structures are not perfectly pure point diffractive, and in order to gain further 
insight into the possible structures of materials, and more generally into the whole concept of 
long-range order, it is necessary to widen the scope of this study to include mixed diffraction 
spectra. In particular, this means that one has to consider structures whose diffraction mea- 
sures contain at least some continuous part. Note that distinct structures may have the same 
diffraction, which makes the corresponding inverse problem difficult (and generally unsolvable 
without further information). In fact, this gets worse outside the realm of pure point spectra, 
and any improvement requires a better understanding of the diffraction of structures with 
some form of randomness. 
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However, when it comes to mixed spectra, relatively little is known, although there are 
many particular examples [23j HI E3 SOI [33 E3 EI] • Even deterministic sets can have mixed 
diffraction spectra, and once any randomness is introduced, this is the norm. Determining 
the exact nature of the diffraction is usually difficult and often simply not known. No doubt, 
the possibilities, both in Nature and in mathematics, for structures with long-range order are 
well beyond what we have presently imagined. This is also made apparent by systems such 
as the pinwheel tiling, compare [49J and references therein, which looks like an amorphous 
structure in diffraction, in spite of being completely regular. In particular, except for the 
trivial Bragg peak at 0, there is no pure point part in the diffraction. 

The primary goal of this paper is to show how the techniques from the theory of stochastic 
point processes may be used when systems with long-range order are subjected to modifica- 
tions in the form of stochastic perturbations. In as much as our primary goal is the study 
of long-range order, the types of point processes of most interest to us are quite different 
from those usually studied in stochastic geometry. For instance, as mentioned above, a phe- 
nomenon of fundamental importance in long-range order is the appearance of Bragg peaks 
in diffraction, which refers to a non-trivial pure point part of the diffraction measure. For 
ergodic point processes, this implies the existence of non-trivial (dynamical) eigenvalues and 
thus excludes weak mixing [62]. Hence, we are particularly interested in systems that lie 
between ergodic and weak mixing. Furthermore, a number of basic and influential mathe- 
matical models of long-range order are deterministic. Even so, the theory of point processes 
is relevant [28] and yields considerable insight. 

To elaborate on this a little, consider the classic Penrose tilings [IT]- Eix a set of two 
generating Penrose prototiles (with matching rule markers) and a set of overall orientations 
for them (so each prototile comes in 10 distinct orientations) . The resulting set of admissible 
tilings of the plane, even if one vertex of the first tile laid down is fixed, is uncountable. 
Replacing each of these tilings by its corresponding vertex point set and allowing all (global) 
translations, collectively we obtain the set X of all Penrose point sets in the given orientation. 
The set X, which is also called the hull, has a natural topology in which it is compact, and it 
carries a unique (and hence ergodic) stationary Borel probability measure \x. The pair (X, \x) 
can now be viewed both as a dynamical system and as a stationary ergodic point process, 
which permits the powerful tools from both subjects to be applied. 

In particular, the diffraction and the dynamical spectrum are linked by the dynamics [6], 
while the diffraction is also the Fourier transform of the first moment measure of the Palm 
measure of the point process |28j . This type of scenario, which applies to many models of 
long-range order, deterministic or otherwise, has not yet been investigated in much detail. It 
would seem desirable then, as a first step, to establish methods, capable of being explicitly 
computable, that would cover typical and much-studied situations and also suggest ways in 
which to generalise what is known, and even move into yet unexplored territory. 

As already implicitly mentioned, mathematical diffraction theory is based on the approach 
set out by Hof in [321 133] . namely via autocorrelations and their measures. Our paper is 
primarily guided by examples, set in as great a generality as we can manage without becoming 
too technical. The examples are selected under the consistent theme that they are computable, 
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while they unify and extend existing results in a systematic way. Briefly, the types of situations 
that we consider are these: 

(i) renewal processes on the real line (Sec. 3), as a versatile, elementary approach to 
one-dimensional phenomena; the main result here is Theorem [TJ 

(ii) randomisation of a given point set A (with a certain discreteness restriction) whose 
diffraction is known, by complex, identically distributed, finite random measures that 
are independently centred at each point of A (Sec. 4); see Theorem [2 

(hi) randomisation of a random point process ^ (with known law) by identically dis- 
tributed, finite, random measures (positive or signed) which are independently centred 
at each point of any realisation of <& (Sec. 5.1-5.4); the main results are formulated in 
Theorems [3] and HI 

(iv) equilibria of critical branching Brownian motions (Sec. 5.5, with Theorem [6]). 

While (i) and (ii) have a bit of a review character, (iii) and (iv) are new in this context. 
The results for (i) and (ii) are included with details and several examples because they have 
immediate applications to practical diffraction analysis, while the results are scattered over 
the literature or only covered implicitly. Since the transfer of methods and results from 
stochastic geometry to mathematical diffraction theory requires a somewhat unusual view 
on measures and their Fourier transforms, we begin with a brief recapitulation of concepts 
needed here. For the sake of completeness and readability (as well as lack of reference), some 
details on the ergodic theorem that we need for (iii) are added as an appendix. 

2. Some recollections from Fourier analysis and diffraction theory 

Throughout the paper, we need various standard notions and results from Fourier analysis 
and measure theory, which we summarise here, introducing our notation at the same time. 

First, let [i be a finite, regular (and possibly complex) Borel measure on M d ; compare |51] 
for background. Its Fourier (or Fourier-Stieltjes) transform is a uniformly continuous function 
on R , defined by 



see |52j for details. This definition includes the Fourier transform of arbitrary Schwartz func- 
tions or integrable functions (the corresponding spaces being denoted by S(R d ) and L 1 ^)) 
by viewing them as Radon-Nikodym densities for Lebesgue measure A, hence as finite mea- 
sures. In this version of the Fourier transform, with the factor 2ir included in the exponent, 
there is no need for prefactors (though the factor 2tt reappears up front under differentiation). 
In particular, one has the usual convolution theorem in the form jT*u = flu, where 



for continuous functions g 6 Co(M d ) (note that we identify finite, complex, regular Borel 
measures on ~K d with continuous linear functionals on Co(M d ), in line with the Riesz-Markov 
representation theorem |51} Thm. 6.19]). 
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Below, we need to go beyond the situation of finite measures. For the introduction of 
unbounded measures, the linear functional point of view is advantageous. Here, an unbounded 
measure is thus understood as a linear functional on the space C c (IR d ) of continuous functions 
with compact support such that, for any compact set K C M d , there is a constant a = a K 
with \^{g)\ < a^lblloo for all continuous g with support in K. The corresponding space 
M.{R. d ) is thus equipped with the vague topology, see [22], Ch. XIII] for details. As before, we 
can identify these measures with the locally finite, complex, regular Borel measures on M. d , 
by an appropriate version of the Riesz-Markov representation theorem; see \12\ Thm. 69.1] 
for a formulation for positive measures and use the polar representation [22, Thm. 13.16.3] 
for an extension to complex measures. The absolute value |^| of fi, also known as the total 
variation measure of /i, is the smallest positive measure such that \n(g)\ < \fi\(\g\) holds for 
all g £ C c (R d ). 

When an unbounded measure fi also defines a tempered distribution, via fJ,((p) = j Rd tpd/j, 
for ip £ 5(IR d ), it is called a tempered measure. Its Fourier transform (as a distribution) is 
then defined via fx{<p) = fi((f>) as usual [50] . so that fi is a tempered distribution. Below, 
we only consider situations where fi is also a measure, hence a linear functional on C c (M. d ). 
Recall that a (complex) measure /U is called translation bounded when 

sup \n\(t + K) < oo 

teR d 

holds for arbitrary compact sets K C ~M. d . Translation boundedness is a sufficient (though 
not a necessary) criterion for a measure to be tempered, see [5U] for details. In this setting, 
we call a measure transformable when it is tempered and when its Fourier transform (as 
a distribution) is again a measure. Transformability of a measure is a difficult question in 
general; see |26j and references therein for details. 

A measure [i G .M(R rf ) is called positive definite, when (i(g*g) > holds for every function 
g E C c (R d ); here, g is the function defined by g(x) = g(—x). Positive definite measures have 
various nice properties, some of which can be summarised as follows; see [13\ Sec. 4] for details 
and proofs. 

Lemma 1. For fi £ _M(M rf ) ; the following properties hold. 

(i) // /i is positive and positive definite, it is translation bounded; 

(ii) // fi is positive definite, it is Fourier transformable, and fi is a positive, translation 
bounded measure; 

(iii) A transformable measure fi is positive definite if and only if fi is a positive measure; 

Moreover, the mapping \i i— > fi defines a bijection between the positive definite and the trans- 
formable positive measures on M. d . □ 

Let us consider some examples that will reappear later. If r C M rf is a lattice (meaning a 
discrete subgroup of ~K d with compact factor group ~K d /r), we write 5p ■= X^er^z f° r the 
corresponding Dirac comb, with 5 X the normalised point measure at x. It is well-known that 
5r is a tempered measure, whose Fourier transform is again a tempered measure. The latter 
is explicitly given by the Poisson summation formula (PSF) in its version for lattice Dirac 
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combs US Ex. 6.22], 

(1) S r = dens(r)5 r *, 

where T* := {x S M rf | x ■ y G Z for all y G T} is the dua/ lattice of -T; see [16J for details. The 
density of T is well-defined and given by dens(i~') = l/|det(i~')|, where det(-T) is the oriented 
volume of a (measurable) fundamental domain of T. It can most easily be calculated as the 
determinant of a lattice basis. Observing |det(P*)| = l/|det(-T)|, a more symmetric version 
of the PSF reads 



(2) (V|det(r)| 6 r ) = V|det(r*)| 6 r *. 

In particular, one has S Z d = 5%d, so that the lattice Dirac comb of Z d is self-dual in this sense. 

Remark 1. Radially symmetric PSF. As an aside of independent interest, let us recall 
the following related formula for a radially symmetric situation in which emerges from 
a simplified model of powder diffraction j3]. Let T and T* be as before, and let r) r {r) and 
T]p*(r) denote the numbers of points of T and T* on centred spheres dB r (0) of radius r. 
The (non-zero) numbers r] r (r) are also called the shelling numbers of the lattice T. If fj, r 
denotes the uniform probability measure on dB r (0), with n = S Q , one has the following 
radial analogue of the PSF in (pQ), 

(3) ( *7r(r) A*r) = dens(r) ^ Vr*( r )Vr, 

rev r ' r£D r , 

where V r = {r > | n r (r) > 0} and analogously for V r ,, see [3] for a proof and further 
details. The formula can also be brought to a more symmetric form, as in Eq. ([2|). 

Another simple, but important, pair of mutual Fourier transforms follows from the relations 
So = A and A = 6q, with A being Lebesgue measure, so that we have 

(4) (<5 + A)" = 5 + A. 

We shall meet this self-dual pair of measures below in Examples [T] and [9j in connection with 
the Poisson process. 

A little less obvious is the following result. 

Lemma 2. Let A denote Lebesgue measure on M. d and < a < d. The function x i— » l/|x| rf-a 
is locally integrable and, when seen as a Radon-Nikodym density for A, defines an absolutely 
continuous and translation bounded measure on M. d . This measure satisfies the identity 

\d-rv W — a 



where the transformed measure is again translation bounded and absolutely continuous. More- 
over, both measures are positive and positive definite. 

Proof. Local integrability of both measures on M rf rests upon that of their densities around 0, 
which follows from rewriting the volume element in polar coordinates, dA(x) = r d ~ 1 drd£2, 
with dJ2 the standard surface element of the unit sphere in M d . Absolute continuity and 
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translation boundedness are then clear, while the Fourier identity follows from a calculation 
with the heat kernel, see [Ml Sec. 2.2.3]. As both measures are clearly positive, they are also 
positive definite by the Bo chner- Schwartz theorem [50\ Thm. IX. 10], compare Lemma [TJ □ 

Incidentally, dividing the identity in Lemma [2] by T(a/2)/ir a / 2 shows that 

a 

( 5 ) d-a TZld^ ► ^0 

r(f)vr— W 

in the vague topology, which follows from the corresponding Fourier transforms of the left 
hand side converging vaguely to A. 

Let us now briefly review the concept of the diffraction measure of a complex measure ui 
as the Fourier transform of the autocorrelation 7 of uj. It motivation comes from the physics 
of diffraction [T7j, while its precise mathematical formulation was pioneered by Hof [32] . 

In general, a complex measure ui need not be transformable, and may thus not be a good 
object for harmonic analysis. In view of Lemma[ll it seems appealing to first attach a positive 
definite measure to uj, which is possible as follows. If uj r denotes the restriction of u to the 
open ball B r of radius r around 0, the natural autocorrelation measure 7 = 7^ is defined as 

6) 7 := km , 

r^oo vol(i3 r j 



provided the limit exists. Here, ju denotes the measure given by jZ(g) = fi(g) for g G C C (K ), 
with g as before. If uj is translation bounded, the one-parameter family of finite measures 
{ to1(*b ' ) I r > 0} is uniformly translation bounded and hence precompact in the vague topol- 
ogy by [32l Prop. 2.2]. One can thus always select converging subsequences to define an 
autocorrelation (which then depends on the sequence of averaging sets). As long as balls 
are used, one speaks of natural autocorrelations. More generally, one may work with any 
averaging sequence 

A = {A n I n G N} 

of relatively compact, open sets A n C M. d that satisfy A n C A n+ \ for all n G N together with 
UneN An = Again, for translation bounded measures u, the corresponding limit in © 
exists, at least along suitable subsequences. 

An important further ingredient is the concept of a van Hove sequence, which is an aver- 
aging sequence with a restricted 'surface to volume' ratio. To formalise this, let K, C C M rf 
be compact and define 

(7) d K C:= ((C + K)\C°) U ((W\C-K)r\C), 

which may be viewed as a ivT-thickened boundary of C. Then, A is called van Hove when, 
for every compact K C M d , 

n— >oo vol[A n ) 

Now, the comparison of limits taken along different averaging sequences makes sense, and 
becomes independent of A for ergodic systems; compare |55j Lemma 1.1]. Also, as follows 
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from [55\ Lemma 1.2], translation bounded measures satisfy the relation 

/n s ,. U n *uT n UJ n *U 

(9) lim — = lim 



oo Vol(^4 n ) n— >oo Vol(A n ) ' 

provided that A is van Hove and one of the limits exists. Here, u n = u\ A , and uj n * uj is 
well-defined by |13^ Prop. 1.13]. This freedom will be used several times below. 

The general situation for a translation bounded measure u is as follows. The van Hove 
property of A implies that \uj\ (A n ) < c vol(A n ) with a constant c > 0. An obvious modification 
of |32} Prop. 2.2] in conjunction with [55\ Lemma 1.2] then gives the following result. 

Lemma 3. Let uj be a translation bounded measure, and A a van Hove averaging sequence. 

With 7n : = * nd 7 n;mo d : = ^p!) ' the famtHeS ^n\n& N} and { 7n;mod | n £ N} are 

uniformly translation bounded and hence precompact in the vague topology. Any accumulation 
point of either family, of which there is at least one, is also an accumulation point of the other 
family, and a translation bounded, positive definite measure. □ 

Lemma [1] applies to any autocorrelation measure, and the corresponding measure 7 is then 
a positive, translation bounded measure. It is called the diffraction measure of uj, relative 
to the averaging sequence A. In ergodic situations, we have no dependence on A and thus 
suppress it. Then, the diffraction measure is also related to the Bartlett spectrum known 
from stochastic geometry, though there are important differences to be discussed later; see 
Remark [15] below. 

In general, an interesting initial question concerns the spectral type of 7, which follows 
from the spectral decomposition 

(10) 7 = (7) pp +(7) sc +(7) ac 

of 7 into its pure point, singular continuous and absolutely continuous parts relative to A, the 
latter being the Haar measure on M d with A([0, l] rf ) = 1. Lattices and regular model sets [551P3] 
are examples with 7 = ( 7 ) pp , while the Thue-Morse and the Rudin-Shapiro sequence show 
singular continuous and absolutely continuous components, respectively; compare [35] and 
references given there. Absolutely continuous components appearing as a result of stochastic 
influence are the main theme below. 

3. Renewal processes in one dimension 

An illustrative class of examples is provided by the classical renewal process on the real 
line, defined by a probability measure g on M + = {x > 0} of finite mean as follows. Starting 
from some initial point, at an arbitrary position, a machine moves to the right with constant 
speed and drops a point on the line with a random waiting time that is distributed according 
to g. When this happens, the clock is reset and the process resumes. In what follows, we 
assume that both the velocity of the machine and the expectation value of g are 1, so that 
we end up (in the limit that we let the initial point move to —00) with realisations that are 
almost surely point sets in R of density 1. 

Clearly, the process just described defines a stationary process. It can thus be analysed 
by considering all realisations which contain the point 0. Moreover, there is a clear (distri- 
butional) symmetry around this point, so that we can determine the autocorrelation (in the 
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sense of ©) of almost all realisations from studying what happens to the right of (we will 
make this approach rigorous in Proposition [T] below). Indeed, if we want to know the fre- 
quency per unit length of the occurrence of two points with distance x (or the corresponding 
density), we need to sum the contributions that x is the first point after 0, the second point, 
the third, and so on. In other words, we almost surely obtain the autocorrelation 

(11) 7 = 5o + v + V 

with v = g + g*g + g*g*g + . . . = YlnLi 6* n an d v as defined above, provided that the sum in 
Eq. (jlip converges properly. Note that the point measure at simply reflects that the almost 
sure density of the resulting point set is 1. In the slightly more general case of a probability 
measure g on M + U {0}, one has the following convergence result. It is essentially a measure 
theoretic reformulation of the main lemma in |25|, Sec. XI.l], but we prefer to give a complete 
proof that is adjusted to our setting. 

Lemma 4. Let g be a probability measure on R + U {0} ; with g(R+) > 0. Then, v n : = 
g + g*g+ • • • + g* n with n £ N defines a sequence of positive measures that converges towards 
a translation bounded measure v in the vague topology. 

Proof. Note that the condition f?(M+) > implies < £>({0}) < lj hence excludes the case 
g = 8q. When g = 5 a for some a > 0, one has v n = J2m=i^ma by a simple convolution 
calculation, and the claim is obvious. In all remaining cases, it is possible to choose some 
a £ M + with g({a}) = and < g([0, a)) = p < 1, so that also g([a, oo)) = 1— p < 1. Since the 
sequence v n is monotonically increasing, the claimed vague convergence follows from showing 
that limsup n „ >00 f n ([0, x)) is bounded by C\ + C2X for some constants Cj. As there are at 
most countably many points y with g({y}) > 0, it is sufficient to show these estimates for all 
x £ M+ with g({x}) =0. In a second step, we then demonstrate that Q* n ([b, b + x)) is 

bounded by 1 + C\ + C2X, independently of b, which establishes translation boundedness. 

If (Xi) iGN denotes a family of i.i.d. random variables, with common distribution according 
to g (and thus values in M + U {0}), one has 

F(X 1 + ... + X m <x) = g* m ([0,x)). 

On the other hand, for the a chosen above, one has the inequality 

[x/a] . v 

F(Xt + ... + X m <x) < P(card{l < i < m \ X, t > a} < x/a) = ^ (7) i 1 ~ pY P m ~^ 



where (™) = whenever £ > m. Observing Y^m=iP m = p/0- ~ P) an d 

e( , ;)(i- P )V-' = (i- P )'~&- 

m=l ^ ' m=0 

for all £ > 1, the previous inequality implies, for arbitrary n £ N, 



l-p 



m=l £=0 v 7 ^ r \ r 1 

which establishes the first claim. 
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For the second estimate, we choose b > 0, x > and observe 

oo oo 

J2s* n [[b, b + x)) =J2 ]P {b<X 1 + ---+X n <b + x) 

n=l n=l 
oo n 

= J2 ^2 F ( X i + ■■■ + x k-i < b < X x + ■ ■ ■ + X k and b < X 1 + • • • + X n < b + x) 

n=l k=l 

oo n 

< Y, ¥ ( Xi + ■■' + ^ b ^ Xi + • • • + Xk ) F ( x k+i + ---+x n<x ) 

n=l k=l 

oo oo 

= J2 F ( X 1 + ■■■ + X k-i < b < Xi + ■ ■ ■ + X k ) ^P(X fc+1 + • • • + X n < x) 

k=l n=k 

oo 

= 1+ ^p(X! + ---+X m <x), 

m=l 

with the convention to treat empty sums of random variables as 0. The last step used the 
i.i.d. property of the random variables together with P(0 < X) = 1 and 

oo 

J>(Xi + • • • + X fc _i <b<X 1 + --- + X k )=l. 

k=l 

In conjunction with our previous estimate, this completes the proof. □ 

When £?({0}) > 0, we are outside the realm of (renewal) point processes, and formula (fTT]) 
for the autocorrelation no longer applies. This case might nevertheless be analysed with the 
methods of Sections 0] and El see Example [7] and Corollary Q] in particular. For the remainder 
of this section, we assume £>({0}) = 0, so that g is a measure on M + ; see Remark [H] below 
for an alternative approach via random counting measures, or |25|, Ch. XI. 9] . 

Proposition 1. Consider a renewal process on the real line, defined by a probability measure 
g on IR + with mean 1. This defines a stationary stochastic process, whose realisations are 
point sets that almost surely possess the autocorrelation measure 7 = 5q + v + u of (jlip . 

Here, v = YlnLi 6* n ^ s a translation bounded positive measure. It satisfies the renewal 
equations 

v = q+q*v and (1 — q)u = g, 

where g is a uniformly continuous function on R. In this setting, the measure 7 is both 
positive and positive definite. 

Proof. The renewal process is a classic stochastic process on the real line which is known to 
be stationary and ergodic; compare [25l Ch. VI. 6] for details. Consequently, the measure 
of occurrence of a pair of points at distance x + dx (or the corresponding density) can be 
calculated by fixing one point at (due to stationarity) and then determining the ensemble 
average for another point at x + dx (due to ergodicity). This is the justification for the 
heuristic reasoning given above, prior to Eq. (jlip . 

By Lemma HJ v is a translation bounded measure, so that the convolution g * v is well 
defined by [13} Prop. 1.13]. The first renewal identity is then clear from the structure of v 
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as a limit, while the second follows by Fourier transform and the convolution theorem. The 
autocorrelation is a positive definite measure by construction, though this is not immediate 
here on the basis of its form as a sum, see [1] for a related discussion. □ 

Let us now consider the spectral type of the resulting diffraction measure for the class of 
point sets generated by a renewal process. This requires a distinction on the basis of the 
support of g. To this end, the second identity of Proposition Q] is helpful, because one has 

< 12 > w - jlgsj 

at all positions k with g(k) 7^ 1. This is in line with summing v as a geometric series, which 
gives the same formula for v(k) for all k with \g(k)\ < 1 and has (|12p as the unique continuous 
extension to all k with \g(k)\ = 1 ^ g(k). In fact, one sees that v(k) is a continuous function 
on the complement of the set {k G R | g(k) = 1}. For most g, the latter set happens to be 
the singleton set {0}. 

In general, a probability measure ^ on R is called lattice-like when its support is a subset 
of a translate of a lattice, see [27] for details. We need a slightly stronger property here, and 
call \x strictly lattice-like (called arithmetic in |25j) when its support is a subset of a lattice. 
So, the difference is that we do not allow any translates here; see [2] for related results. 

Lemma 5. If fj, is a probability measure on R, its Fourier transform, fi(k), is a uniformly 
continuous and positive definite function on R, with < ju(0) = 1. 

Moreover, the following three properties are equivalent. 

(i) card{/c G R | /x(fc) = 1} > 1; 

(ii) card{£; G R \ fi(k) = 1} = oo; 

(iii) supp(//) is strictly lattice like. 

Proof. One has fi(k) = f R e~ 2nlkx dfi(x) , whence the first claims are standard consequences 
of Fourier analysis; compare [48j Prop. 5.2.1] and [52l Sec. 1.3.3]. 

If fj, = XLer p( x )3x f° r a lattice fcR, with p(x) > and Ylx&r p( x ) = 1> one nas 

m = ^2p(x)e- 2mkx , 

so that fi(k) = 1 for any k G r* . In particular, r* C {k G R | //(A;) = 1}, so that we have the 
implications (iii) (ii) (i). 

Conversely, if fi(k) = 1 for some k ^ 0, one has L e - 27rjfca; d/u,(x) = 1 and hence 

(13) / (1 - cos(2vr/cx)) dfj,(x) = I (l - cos(2ttA;j;)) d(i(x) = 0, 

Jl ^supp(/i) 

where supp(/^), the support of the probability measure [a, is a closed subset of R and mea- 
surable. The integrand is a continuous non-negative function that, due to k ^ 0, vanishes 
precisely on the set ^Z, which is a lattice. 

Write supp(/x) = AUB as a disjoint union of measurable sets, with A = supp(/x) n and 
5 = supp(^)n(R\iZ). We can now split the second integral in (|13p into an integral over A, 
which vanishes because the integrand does, and one over the set B, which would give a positive 
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contribution by standard arguments, unless B = 0. But this means supp(/z) = A C rZ, so 
that (i) =>• (iii), which establishes the result. □ 

At this point, we can state the main result of this section, the diffraction properties of 
renewal processes; compare [191 Ex. 8.2(b)] for a special case. 

Theorem 1. Let g be a probability measure on M + with mean 1, and assume that g is not 
strictly lattice-like. Assume further that a moment of g of order 1 + e exists for some e > 0. 
Then, the point sets obtained from the stationary renewal process based on g almost surely 
have a diffraction measure of the form 7 = 5q + (7) ac with 

!_ |-(Ml2 



where h is a continuous function on E \ {0} that is locally integrable. It is given by 

_ 2 (\g(k)\ 2 ~ Md(k))) 
h[k) ~ \l~Q{k)\ 2 
and measures the difference from a constant background as described by A. 

When g is strictly lattice-like, the pure point part becomes a lattice Dirac comb, and the 
behaviour of h at repeats at each point of the underlying lattice, see Remark [3] for details. 

Proof. The process has a well-defined autocorrelation 7, by an application of Proposition [H in 
the sense that almost every realisation of the process is a point set A with this autocorrelation. 
Since 7 is a positive definite measure, it is Fourier transformable by LemmaHUh), with 7 being 
a positive measure on R. 

The point measure at with intensity 1 reflects the fact that the resulting point set A 
almost surely has density 1. To see this, define g n = and h n = g n * g n . Here, h n is 

a positive definite, tent-shaped function with support [— n, n] and maximal value ^ at 0. It 
has (inverse) Fourier transform h n (k) = ( ""^T^ ) 2 > wmcn is a non-negative function (with 
maximum value 1 at k = 0) that concentrates around as n —* 00. Let oj r := S An , rr y Using 

j6|) together with (w r * (ly ) (g n * g~ n ) > 0, it is not difficult to see that 7(/i n ) n ^°° > (dens(yl)) 2 , 
which is almost surely 1 (for this, assume first that r > n > 1, then take the limit r — > 00 
followed by the limit n — > 00). On the other hand, one has 

1{K) =%hn) =7(M 7(W), 

due to the concentration property of h n (in particular, for all e > 0, one verifies the relation 
7(-B e (0)) > (dens(yl)) > 0, which proves the existence of a point measure at 0). 

Due to the assumption that supp(^) is not contained in a lattice, we may invoke Lemma [5] 
to see that g(k) 7^ 1 whenever k 7^ 0, so that we have pointwise convergence 

onM \ {0}, and similarly for v. Since g is uniformly continuous on M and g(k) / lonK\ {0}, 
both v and v are represented, onR\ {0}, by continuous Radon-Nikodym densities. Writing 
(5q + v + v ) = (1 — h)X, hence {y + u) = —h\, the formula for h follows from v = v. 
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It remains to show that 1 — h is locally integrable near 0. Let X be a random variable with 
distribution g. Since the latter has mean 1 and our assumption guarantees that (X 1+£ ) = 
J °° x l+s dg(x) < oo, we have the Taylor series expansion 

g(k) = 1 -2irik + 0(\k\ l+£ ), as \k\ -► 0, 

by an application of |614 Thm. 1.5.4]. Inserting this into the expression for h results in 

h(k) = 2 + <D(k~ 1+£ ), as \k\ -► 0, 

which establishes integrability around 0, and thus absolute continuity of the measure (1 — h)X. 

As the contribution to the central peak is already completely accounted for by the term 
So, the claim follows. □ 

Remark 2. Asymptotic behaviour of h. When, in the setting of Theorem [H the second 
moment of g exists, one obtains from |6H Thm. 1.5.3] the slightly stronger expansion 

g(k) = 1 - 2nik - 2vr 2 (A 2 ) k 2 + o(\k\ 2 ) , as \k\ -> 0. 

This leads to the asymptotic behaviour 

h(k) = 2- (X 2 ) +0(1), as|A;|^0, 

which implies that h is bounded and can continuously be extended to h(0) = 2— (X 2 ) = I — a 2 , 
where a 2 is the variance of g. Clearly, the existence of higher moments implies stronger 
smoothness properties. 

Remark 3. Complement of Theorem CD When g happens to be strictly lattice like, the 
Z-span of the finite or uniformly discret^l set supp(^) is a lattice of the form r = 6Z, where 
b > is unique (in other words, r is the coarsest lattice in M that contains supp(^)). Then, 
one finds the diffraction 

7 = h/b + (1 " h)X, 

with the function h from Theorem [TJ Note that h is well-defined (and continuous) on R \ r*, 
with r* = TLjb being the dual lattice of I '. Moreover, it is locally integrable around all 
points of r*, so that (1 — h)X is again an absolutely continuous measure. Note that, since 
the underlying point set is always a subset of r, the diffraction measure is periodic, with r* 
as its lattice of periods; compare [2] for general results in this direction. 

When supp(g) is a finite set, one is in the situation of a random tiling with finitely many 
prototiles. A more detailed discussion, together with an explicit calculation of h for this case, 
is given in [U Thm. 2]; see also Example [5] and Remark [71 

Let us turn to some examples, for which we employ the Heavyside function, 

1, if x > 0, 

(11) 9(x) := { i, ifx = 0, 

0, if x < 0. 



^Recall that a set S C R d is called uniformly discrete when there is a number s > such that the distance 
between any two distinct points of S is at least s. 
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This formulation of © has some advantage for formal calculations around generalised functions 
and their Fourier transforms. 

Example 1. PoiSSON PROCESS ON the real line. The probably best-known stochastic 
process is the classical (homogeneous) Poisson process on the line, with intensity 1, where 
q = f\ is given by the density 

f(x) = e- x 0(x). 

It is easy to check that the convolution of n + 1 copies of this function yields e~ x x n 0(x)/n\, 
which gives v = OX. As the intensity is 1, this results in the autocorrelation 

7 = 5q + u + T> = 5 + A 

and thus in the diffraction 7 = 7, compare Eq. 

Remark 4. Characterisation of Poisson processes. Let N denote a homogeneous 
Poisson process on the real line, so that, for any measurable A C M, N(A) is the number of 
renewal points that fall into A. It is well-known that N(A) is then Poisson-(A(^4))-distributed, 
which means that 

e -KA) (\(A)) k 

F(N(A) = k) = { { }) 

k\ 

with k E No, and that, for any collection of pairwise disjoint sets A\, A2, ■ ■ ■ , A m , the random 
numbers N(A{), N(A m ) are independent. In fact, this property characterises the Poisson 
process (compare |19^ Ch. 2.1]), and it can serve as a definition in higher dimensions or in 
more general measure spaces, to which the renewal process cannot be extended. 

Example 2. Renewal process with repulsion. A perhaps more interesting example in 
this spirit is given by the density 

f{x) = Axe- 2x 0{x). 

It is normalised and has mean 1, as in Example [H but models a repulsion of points for small 
distances. Note that this distribution can be realised out of Example [1] by taking only every 
second point, followed by a rescaling of time. 

By induction (or by using well-known properties of the gamma distributions, compare 
Sec. II.2]), one checks that 

= (2^171 x 2n ^e- 2x e(x% 
which finally results in the autocorrelation 

7 = J + (i- e - 4 W)A = 5 + A - e -4 !*! A 
and in the diffraction measure 

2 + (irk) 2 x 2A 



4+ (irk) 2 u 4+(vrA:) 2 ' 

This is illustrated in Figure [TJ The 'dip' in the absolutely continuous part around 0, and 
thus the deviation from the previous example, reflects the effectively repulsive nature of the 
stochastic process when viewed from the perspective of neighbouring points. 
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Figure 1 . Absolutely continuous part of the diffraction measure from Exam- 
ple EJ for a = 0.7 (top curve), a = 1 (horizontal line, which also represents 
Example CQ , a = 2 (see also Example [2]) and a = 8 (overshooting curve) . 



Example 3. Renewal process with gamma law of mean 1. The previous two examples 
are special cases of the gamma family of measures. For fixed mean 1, they are parametrised 
by a real number a > via g Q = f a X and the density 



(15) 



fa(x) 



a 



r(a) 



x a-l e -ax 



0(x) 



While a = 1 is the 'interaction-free' Poisson process of Example [H the density implies an 
effectively attractive (repulsive) nature of the process for < a < 1 (for a > 1). When 
a = k £ N, the process can also be interpreted as a modified Poisson process where one keeps 
only every kth point (followed by an appropriate rescaling). 

Observing /^ n (x) = r "^) x na_1 e~ ax 0{x) for n G N, this leads to the measure 



(16) 



g a 0X 



with 



9 a {x) 



a e 



y> (ax)""' 1 
^ T(na) 



n=l 



9a( x ) = 1- The calculations result in the autocorre- 



Note that, for fixed a, one has lim^^oo 
lation 

la = S + g a {\x\) A 

and in the diffraction 7 Q = 5q + (1 — h a ) A, where h a is the symmetric function defined by 

= 2(l-Re((l + 2^fc/«)°)) 
|l - (1 + 2mk/a) a \ 2 

The latter follows from the general form of h in Theorem [TJ together with the observation 
that f a (k) = (1 + 2mk/a)~ a . 

It is easy to see that lim^ — h a (^k^) — 0, for any fixed ol ]> 0, which makes the role of h a as 
the deviation from the Poisson process diffraction more transparent, where a = 1 and hi = 0. 
Note also that lima^oo = 5i in the vague topology, in line with the limits mentioned 
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before. This can nicely be studied in a series of plots of the diffraction with growing value of 
the parameter a. Figure [1] shows some initial cases. 

Remark 5. Construction of Delone sets. Of particular interest in the applications 
are Delone sets (which are point sets that are both uniformly discrete and relatively dense), 
because points (representing atoms, say) should neither be too close nor too far apart. Such 
sets can also arise from a renewal process. In fact, if one considers a probability measure g on 
M+, the resulting point sets are always Delone sets when supp(^) C [a, b] with < a < b < oo, 
and conversely. This equivalence does not depend on the nature of g on [a, b], while the local 
complexity of the resulting point sets does. In particular, if g is absolutely continuous, the 
point sets will not have finite local complexity (see below for a definition). 

It is clear that no absolutely continuous g is lattice-like, hence certainly not strictly lattice- 
like, so that all these examples match Theorem [TJ Probability measures g with supp(^) 
contained in a lattice are covered by Remark [3l They are of interest because they form a link 
to point sets and tilings of finite local complexity, which have only finitely many patches of a 
given size (up to translations). Let us consider some examples. 

Example 4. Deterministic lattice case. The simplest case is g = 8\. From 5i*5\ =62, 
one sees that v = 5 N and hence 

7 = ^0 + <5n + <5-n = h > 
which is a lattice Dirac comb, with Fourier transform 

7 = h 

according to the Poisson summation formula ([1]) . This is the deterministic case of the integer 
lattice Z, covered in this setting. 

Remark 6. Deterministic limit of Example [3j The last example can also be seen as a 
limiting case of the measure g a defined by Eq. (fT5l) . In particular, one has linia^oo g a = 5\ 
and liuLo^oo v a = 5j$, with v a as in (|16p and both limits to be understood in the vague 
topology. This can also be seen by means of the strong law of large numbers. For each 
n E N, by well-known divisibility properties of the family of Gamma distributions, g n is the 
distribution of 



1 

/ . X i 
i=l 



where the Xi are independent and exponentially distributed random variables with mean 1. 
This sum then concentrates around 1, with a standard deviation of order 1/y/n. 

Example 5. Random tilings with finitely many prototiles. Consider the measure 

g = a5 a + (l-a)8b, 

with a E (0, 1) and a, b > 0, subject to the restriction aa + (l — a)b = 1 to ensure density 1. 
Each realisation of the corresponding renewal process results in a point set that can also be 
viewed as a random tiling on the line with two prototiles, of lengths a and b. As before, place 
a normalised point measure at each point of the realisation. Then, the diffraction (almost 
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surely) has a pure point and an absolutely continuous part, but no singular continuous one. 
The pure point part can be just 5o (when b/a is irrational) or a lattice comb (see Remark [3]); 
details are given in [JJ, including an explicit formula for the AC part. 

This has a straight- forward generalisation to any finite number of prototiles, with a similar 
result. Also in this case, there is an explicit formula for the diffraction measure, which was 
derived in [3j by a direct method, without using the renewal process. 

Remark 7. Continuous diffraction with 'needles'. Looking back at Lemma El one 
realises that Example [5] revolves around the lattice condition in an interesting way. Namely, 
even if q is not strictly lattice-like, supp(p) for a random tiling example with finitely many 
prototiles is a finite set, and thus a subset of a Meyer set (which is a relatively dense set A 
whose difference set A— A is uniformly discrete) . We then know from the harmonic analysis of 
Meyer sets, compare |44| and references therein, that g(k) will come e-close to 1 with bounded 
gaps in k. This means that the diffraction measure, though it is absolutely continuous apart 
from the central peak at k = 0, will develop sharp 'needles' that are close to point measures in 
the vague topology — a phenomenon that was also observed in [4] on the basis of the explicit 
solution. 

4. Arbitrary dimensions: Elementary approach 

Let us now develop some intuition for the influence of randomness on the diffraction of point 
sets and certain structures derived from them in Euclidean spaces of arbitrary dimension. In 
this section, our point of view is from a single point set A C M rf that is being modified 
randomly, by replacing each point by a complex, finite, random cluster. This is still relatively 
easy as long as A is sufficiently 'nice'. In Section [5j we revisit this situation from the point of 
view of a stationary ergodic point process, which treats almost all of its realisations at once 
and permits a larger generality for the sets A, though the clusters will then be restricted to 
positive or signed measures. 

Let A C M rf be a fixed point set, which we assume to be of finite local complexity (FLC). 
By definition, this means that there are only finitely many distinct patches of any given 
size (up to translations) in A. This property is equivalent to the difference set A — A being 
locally finite [55], the latter saying that K D (A — A) is a finite set for all compact K C M d . 
In particular, since is then isolated in A — A, the set A itself is uniformly discrete; see 
Remark 1111 for a possible extension. Attached to A is its Dirac comb 5 a = X^e/l ^ x > w hich is 
a translation bounded measure, as a consequence of the FLC property. We associate to 5a the 
autocorrelation and the diffraction measure as explained in Section for a suitably chosen 
averaging sequence A = {A n \ n G N} of van Hove type. A natural choice is A n = B rn (0), 
with B r (0) denoting the open ball of radius r around 0, for a non-decreasing series of radii 
with r n "^°°> oo (alternatively, nested cubes are also quite common). 

Set A n = An A n (so that A n SA'm the obvious local topology [55]) and consider 
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We now make the assumption that the limit 

(17) lim j A =: j A 

n— ->oo ' 

exists in the vague topology, which is then the autocorrelation measure of the set A relative 
to the averaging sequence A. 

Remark 8. Accumulation points versus limits. Due to translation boundedness of 
S A , the sequence of measures j A n always has points of accumulation; see (32^ Prop. 2.2] and 
Lemma [3l Consequently, one can always select a subsequence of A for which the assumption 
(|17p is satisfied. This remains true even if we relax the nesting condition for A. In this sense, 
when the autocorrelation is not unique (as in the example of the visible lattice points without 
nesting [9J), we simply select one of the possible autocorrelations by a suitable choice of A. 
Our results below apply to any autocorrelation of this kind separately. In this sense, the 
assumption made in (fT7"|) is not restrictive. 

As briefly explained in Section [21 see Lemma [3l the van Hove property of A in the context 
of (|17p implies that one also has 

(18) lim "f An . mod := lim ^"*/^ = 1 A , 



the difference between the two approximating measures in (|17p and (118p being a 'surface term' 
that vanishes in the infinite volume limit n — > oo. Eq. (|17|) explicitly shows that the measure 
j A is positive definite (hence transformable by Lemma [1]), while ()18p is easier to work with 
for (pointwise) calculations in the presence of random modifications as introduced below. 

Since A — A is locally finite by assumption, Eq. (|18p is equivalent to the existence of all the 
pointwise limits 

(19) lim rj n (z) =: r/(z), 



with the approximating coefficients 

Vn( z ) z 



card{x € A n \ x — z £ A} 



vol(A„) 

where i](z) = for any z A — A. Clearly, the measure ^ A as well as the coefficients rj(z) 
may (and generally will) depend on the averaging sequence; compare Remark [HJ 

The next step consists in modifying A by a random process in a local way. To come to a 
reasonably general formulation that includes several notions of randomness known from lattice 
theory, compare [60\ 163] . we employ a formulation with finite, random, complex measures. 
Let Q denote a measure-valued random variable, and Q the corresponding law, which is 
itself a probability measure on A^bd = -Mbd(^ d )) the space of finite complex measures on 
M. d . To keep the notation compact, we use the symbol Eg for the various expectation values 
that arise in connection with (f2,Q). In particular, we write Eg(l?) = Jj^ bd uj dQ(iv), where 
uj refers to the realisations of Q as usual. Note that we also refer via the index Q to the 
underlying law for one random variable for more complicated expectation values, rather than 
using the underlying (though hidden) probability space. This will be explained in more detail 
in Section [5] below. 
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To proceed, we need a version of the strong law of large numbers (SLLN) for measures. 

Lemma 6. Let (fti)i^ be a sequence of integrable, finite, i.i.d. random measures, with com- 
mon law Q. Then, with probability 1, one has 

1 n 

n t—f 

i=i 

in the vague topology. 

Proof. By definition, integrability means that Eq(|1?i|), which is independent of i £ N, is a 
finite measure. As the space of continuous functions C c (R d ) is separable, the almost sure 
convergence of the measures follows from the almost sure convergence of — ^2^=1 ^(^O f° r 
an arbitrary (but fixed) bounded, continuous function tp. This, in turn, follows from the 
conventional SLLN [24J, possibly after splitting the sums into their real and imaginary parts 
and applying the SLLN twice. □ 

Recall that ui is the measure defined by u) (99) = u)((p). Let ft and ft' be two independent 
random measures, with the same law Q, and such that Eq(|J?|) is a finite measure, and also 
assume the second moment condition Eg((|l7|(IR <:i )) 2 ) < 00. Then, the convolution ft*ft' is 
well defined, and one obtains from elementary calculations the important relations 

(20) E Q (ft) = E^(ft) and E Q (ft*ft') = E Q (ft) * E^ft'), 

the second due to the assumed independence. 

Let us now fix an FLC set A, which is assumed to possess the autocorrelation measure 
j A relative to the van Hove averaging sequence A chosen, and consider the family (ft x ) x£j \ 
of integrable, complex, i.i.d. random measures, with common law Q and subject to the mo- 
ment conditions mentioned above. When ft is any representative of these random measures, 
Eq(|J?|) is a finite measure by assumption, and the measure-valued expectations Eq(]7) and 
Eq(/2 *ft) exist (note that also Eq(\ft*ft\) is a finite measure, due to the condition on the 
second moment). We are now interested in the random object 

(21) =Y,^x*S x , 

which is almost surely a locally finite measure (though not necessarily translation bounded). 

To see this, we observe that, for any bounded Borel set B C M d , the sum J2 xG/ [\(ftx*o~x)(B)\ 
converges almost surely, since 

E Q (\(ft x *6 x )(B)\)=E Q (\ft(B-x)\) < {E Q (\ft\))(B-x) 

and the convolution J yl *Eg(|i?|) is a well-defined locally finite measure due to the translation 
boundedness of S A (note that the summands in ^2 x£y [\ft x * S X \(B) are non-negative, hence 
convergence of the means implies almost sure convergence). As the Borel <7-algebra on R rf 
is countably generated, we can find a set of Q-measure 1 on which the sum (I2ip converges 
(absolutely) for each Borel set B, and the limit is a measure. 
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Let A be fixed and assume for simplicity that each A n is invariant under x >— ► —x. The 

A 



(n-th) approximating autocorrelation of 6^*' reads 



< 22 > « - ^cj <?V - ^ * «V 5>» * M k ■ 

For certain pointwise calculations and arguments, it will be more convenient below to consider 
the modified approximating autocorrelation 

( 23 ) C;mod ^= ^4"t(E ^*^)*(E^- 

To this end, we need a probabilistic analogue of Eq. Q. 

Proposition 2. Almost surely, T^n 0/ ([22]) and J^n-mod 0/ (f23j) define sequences of locally 
finite random measures. Moreover, we can choose a strictly increasing subsequence (nk) keN 
such that, in the vague topology, we almost surely have 

M,Tii. M,n fe ;mod 

In particular, if 7^ or 7^ 

■mod a l m °st surely converges to 7 along A or along a subsequence 
of it, we can choose a subsequence A' of A so that both sequences almost surely converge to 7 
along A' . 

Proof. We abbreviate //(•) := E(|ft|(-)), := E(\f2\(R d )\f2\(-)). Due to the assumption 
E((|i7|(M d )) 2 ) < 00, both \i and v are finite, positive measures. Consequently, since A is 
uniformly discrete, n* 5a and 1/ * 5a are translation bounded by [T31 Prop. 1.13] (and thus 
certainly locally finite). 

Let us first verify that the expression in f|23j) almost surely defines a locally finite measure 
(the estimate for f]22H is completely analogous, with the same upper bound). If S C M rf is a 
bounded Borel set, we have 



vol(^ n )IEQ(l7JS ;m odl( 5 )) < E E 4/ Is^ + ^dd^h^WdOO,!*^)^ 

< E / l B (« + ^)d(M*4)(^)d(ju*<L. 2/ )(v) 

2/6A\{a;} 

+ VEq / l B (u + u)d|n a |(iOd|J&(t;) 

_^ 4 L/K d xR d 



xeA n 



< (i»* 5 A n ) * + ^(M d ) card(A„) < 00. 



Thus, by arguments analogous to those used before, the sum on the right hand side of (|23[) 
almost surely converges absolutely when applied to any bounded Borel set B. Again, since 
the Borel cr-algebra on R d is countably generated, this suffices to ensure that 7^. mod is a 
locally finite, random measure. 
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Furthermore, again for a bounded Borel set B, one has 



(24) 



x,y&A 
x^y 



vol(A0E Q (l7S-7S ;mod l(S)) 

<J2 e q\[ iB^+v)\i A ^)UM-i A jx)\d{\n x \*5 x )(u)d{\n y \*5. y )(v) 

+ IW^W^) - lA n (x)\d[\n x \*6 x )(u)d(\fr x \*6- x )(v) 



xgA 



We estimate the term in the third line of (1241) as follows: 



x,y£A 
x^y 



V / l B (u + v) \l An (u)l An (v) -l An (x)\d(fj,*5 x )(u)d(fi*5- y )(v) 

= EE/ 1 B(« + «)|lA B («)l^(«)-l|d(^*^)(«)d(M**- v )(«) 
+ V V / l B (u + u)l^(«)l x (v)d(/x *$*)(«) <!(/!*$_,,)(«) 



xgA\A n y&A' 
y^x 



< 



+ 



l B (u + v) \l An (u)l An (v) - l\d(fi*5 An )(u)d((x*5 A )(v) 
/ 1b (« + v ) 1a„ («)1a„ («) d(^x * t^J («) d(// * 6 A ) (v) 



< 



where, in the first inequality, we have used the fact that removing the restriction x ^ y in the 
summation only adds positive terms (note that \x is a positive measure by definition), and 
employed the estimate \ l An (u)l An (v) — l| < l K d^ (u) + l R d\ A i v ) for the second inequality. 
There is a constant c B that depends on B (as well as on \x and the averaging sequence) and 
a sequence d n — * that is independent of B such that the sum in the last two lines above is 
bounded by d n c B vol(A n ). This comes from the fact that there are only contributions from 
'surface terms'; compare the arguments in [55j . 



DIFFRACTION OF STOCHASTIC POINT SETS 

By way of example, we verify that, for any R > 0, 



21 



<(m d )\{ 



xeA n \ d(x, R d \A n ) < R}\ + fi(R d \B R )\A n \) sup [n * 6 A ) (B + x) 

' x€R d 



Note that this (together with analogous statements for the two other summands above) yields 
the claim by the van Hove property of the averaging sequence A. 

Observe next that, for a (possibly) complex measure £ with |f | (R d ) < oo, a translation 
bounded measure v and a bounded Borel set B, we have (by an application of [32} Prop. 2.2] 
and its proof) the estimate |(f * v){B)\ < |f \(R d ) sup xeR d\iy\(B + x) < oo. Finally, note that 

H* 5 A is translation bounded by [13, Prop. 1.13], and that 

( M * 5 An ) | (R d ) = £>| * 5 x )(R d \A n )+ Yl (Ml * ^)(M d \A n ) 

</i(M d )|{x e A, I d(s,R d \4,) < +v(R d \B R )\A n \. 



Similarly, the term in the fourth line of (|24f) is bounded from above by 
VEq / l B (u + v) (l Rd \ A n (u) + l md \ A («)) d(|J2 x | *<5 s )(u)d(|f2 a .| *<y_ a ,)(« 
+ Y E Q [ l B (u + v)l A Ju)l A Jv)d{\n x \*5 x )(u)d{\n x \*5„ x )(v) 

x & A\A n U® d x® d 

< ^E Q (|^|(iR d ) / WMddn,!*^)^)) 

+ Y E Q (|^|(M d ) / l RdV4 >)d(|^|*,5_,) 

+ V E Q (|^|(R d ) / l An ( U )d(|^|*4)( U 
*eA\A n V ^ 

(25) = {v*5 A J(R d \A n ) + (i^ n )(R d \A n ) + {i;*5 AXAn )(A n ) < <vol(A n ), 
with a sequence d' n — > 0. Combining the above estimates, we obtain 



< 26) -V vol(A, 

for a sequence d'^ — > 0; hence, for e > 0, 



< (c B + l)< 
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by Markov's inequality. If we choose (nfc) fcgN such that ^ fc d'^ k < oo, we obtain from ([2 
and the (first) Borel-Cantelli lemma that 

(28) rj 2 -. — — > U, almost surely as k — ► oo . 

vol (A* J 

By (f27|) . we may choose the subsequence (nfc) fcgN independently of B in such a way that, for 
each bounded Borel set B, ()28[) holds almost surely. Finally, since the Borel cr-algebra on R rf 
is countably generated, this implies the main claim of the lemma. The last statement is then 
obvious. □ 

Let us now resume our study of j}^ . Invoking Proposition [2] and replacing the averaging 
sequence A = (A n ) n£N by the subsequence chosen there, we may use the modified measures 
7/ln-mod as our approximating measures. Observing 



the modified autocorrelation approximant reads 

-yfflmod = (E * (E^y* 6 - 

( 29 ) /I ~ X 

= E biixr E ^*^-*)*** =: E 

where we now need to analyse the behaviour of the random measures Cz^n- 
Let us first look at z = 0, where we obtain 

(30) C - "ffff J^, E — ^-^(fl.g) (,,) 

by an application of Lemma [6l Note that dens(yl) = r/(0) as introduced in Eq. (|19p . Next, 
assume z *E A — A with z / 0. Then, we split into two sums, 

(31) cii? = ( E (0) ^ * 3— + E (1) ^ * <v*) > 

x—zGA y—zGA 

where the upper index stands for the following additional restriction: Given z, our point set 
A is the disjoint union of countably many maximal linear chains of the form 

{. . . , x + 2z, x + z, x, x — z, x — 2z, . . .} 

with all points lying in A, and x being chosen as its representative. Such a chain may be 
finite or infinite, but has no gaps by construction. For each of these chains, the random 
measures f2 x+mz * (2 x +( m -i)z) are identically distributed, but not independent (due to the 
index overlap). However, those with m even (type (0)) are mutually independent, as are 
those with m odd (type (1)). This way, each element of A inherits the type as a label, and 
the terms in (j3Tj) are distributed to the two sums according to their type. This approach 
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guarantees that the terms for C^n which already showed up in C* n-i en d up in the same sum 
as before, no matter what the detailed structure of the (nested) averaging sequence might be. 
We also split the number of terms 

card{3; G A n \ x - z G A} = + 

accordingly. We can now rewrite our previous expression in the form 

r ^ r{ n) _ cardjx eA n \x-z£A} ( £< > , E W 



V oi(A n ) v 40) + N m N (o) N (o) + N (o) N m 

where the term in brackets is a convex combination of two random measures 

E (0) /M 0) and 

Y^/Njp. By ([19]), the factor in front of the bracket converges to rj(z). When this limit is 
non-zero, we know that n ~~ > °°> 00 for i G {0, 1}, so that Lemma [6] and Eq. (|20[) imply 

(33) ~^lf Eg(fi)»E^(fi) (a.s.). 

Although we do not know whether the rational prefactors in (|32p converge, we have a convex 
combination of two sequences that each almost surely converge to the same limit, which must 
then also be the limit of their convex combination. Put together, this gives 

(34) Cg } i ^ ^)-Eq(«)*e5«) (a.s.) 

for all £ G A - A with z ^ 0. 

These considerations will be sufficient when the random measures almost surely have a 
(deterministic) compact support. To formulate the main result of this section in greater 
generality, we need one further technical property. For brevity, we write = M d \ B r (0). 

Lemma 7. If M' is a set of uniformly translation bounded, positive measures and v a finite, 
positive measure on ~R d , there is a sequence Rk f 00 such that 

00 

V sup (u\ * v){K) < 00 
holds for any compact set K C R . 

Proof. Since v is a finite, positive measure and v{Bf) a decreasing function that tends to 
as r — > 00, we can choose radii with v(B c Rk ) < 1/k 2 , so that Ylfc=l u i^R ) < tt 2 /6 < 00. 
Moreover, we may do this in such a way that the differences between consecutive radii do not 
decrease, meaning that R2 > 2Ri and Rk+2 — Rk+i > Rk+i — Rk for all k G N. 

Uniform translation boundedness of Ai' means that, for any compact K C R , there is a 
positive constant a K with uj{x + X) < a^, simultaneously for all x G M rf and all a; G M' . If 
i<C C M d is compact, we have K C B r (0) for some r > 0. If R > r, one has 



)(#)=/ M* + y)iBc(s)<M*)<MiO 

JR d xlR d 

= / u{(K-y)nB c R )dv(y) < a K u{B R _ 

JR d 
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where the last step follows because ui((K — y) n Bf^ vanishes whenever \y\ < R — r. For radii 
< R < r, the bound is simply given by a K i>(M> d ), which is finite. 

Now, for some m G N, we have R m — r > R\, with the sequence of radii chosen before. The 
additional difference property of the radii makes sure that the radii R m +k with k G N give 
a summable contribution, while the remaining terms are finite by construction. Since this 
argument is uniform in u G M.' and holds for all compact K C M d , our claim follows. □ 

Theorem 2. Let A C M. d be an FLC point set such that its Dirac comb 8 a possesses the 
autocorrelation measure 7^ of (|T~T|) . relative to the fixed averaging sequence A, and thus the 
diffraction measure 7^. Let (f2 x ) xG A be a family of integrable, complex, i.i.d. random measures 
with common law Q and finite second moment measure, with Q being any representative of 
this family, and consider the random measure 8 A °f (EH) • 

Then, possibly after replacing A by a suitable subsequence A' , the sequence of approximating 
measures 7^] of (|22l) almost surely converges, as n — > oo ; to the positive definite, translation 
bounded autocorrelation measure 

1A,Q = (E Q (tf)*E^2)) * lA + dens(A)(E Q (n*f2)-E Q (n)*E^(72))* 8 . 

This measure has the Fourier transform 

1a,Q = \E^)\ 2 -% + den S (A)(E Q (n*fi)-E Q (n)*E^(n))~-\, 

which is the almost sure diffraction measure of the random measure 8]P relative to A' . 

Proof. The previous calculations establish the individual almost sure convergence of the 
(countably many) measures , with the limits as given in Eqs. ([30]) and ([33| . Our as- 
sumptions on Q ensure that Eq{Q) * Eq([2) in (|34j) is a finite, positive definite measure, 
which is concentrated at in the sense that \Eq({2) * E Q (f2)\(R d \ B r (0)) 0, while 

{E Q {Q) * E^[Q)){g * g) > for all ^ g G C c (R d ). In view of (JH and dMD, our (claimed) 
almost sure limit "f^Q inherits translation boundedness from 7^. 

The (deterministic) measure 7^ q is positive definite, hence transformable by LemmaHJ Its 
Fourier transform has the form claimed as a result of the convolution theorem |13] Exes. 4.18]. 
The latter is applicable here because all expectation measures involved are finite measures, 
so that their Fourier transforms are represented by uniformly continuous functions on M. d . 

It remains to establish the limit property. Let us first assume that there is a (deterministic) 
compact set C so that supp(i? x ) C C almost surely. This implies supp(cl^) C C — C for 
all n and z, so that only terms from finitely many z G A — A contribute to 74 „. mo( j on any 
compact K C In this case, we may use an elementary pointwise calculation to see that 
T/ln^-mod tends to the claimed limit, and Proposition [2] gives the assertion. 

In the general case, this simple argument is not conclusive, and we need to estimate putative 
contributions from distant points z G A - A to 7^ 

■ mo d' ^ or bounded K,B B 0, with 
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/*(■) :=E (|^|(0) as above, we have 
(35) 



E 



Q 



y — 

z&A-A v ' x£A n 
z4B x—zGA 



{C2 x *n x - Z *5 Z ){K) 



< 



, „ vol(A t , 

zeA-A y ' xeA n 

z4B x~z£A 



vol(An) E (m^ z )(k) 

z£A-A K n> xeA n 
z4B x—z&A 



1 



Vol^n) 



E p 



zeA-A,x€A n 
x-zeA, z4_B 



A t *M*(7/l,n;mod| 



l d \B' 



(K) =: 4> n {K,B). 



Since {7^ in ;mod I € N} are uniformly translation bounded by Lemma [3J we can choose a 
sequence of radii f oo according to Lemma [7] such that, for all compact K C M. d , 



(36) 



V sup </) n (K,B Rk (0)) < oo. 

nGN 



fc=l 



We have proved above that, for each z G A — A, Q^n ra ^°° > Cz!oo, almost surely in the 
vague topology, where = dens(yl) -Eq(,J?* i?) and Czfoo = ??(<£) -Eq(J7) *Eq(I2) for z / 0. 
Possibly after passing to another subsequence, we may now assume that the convergence is 
so fast that, for any g G C c (K d ), 



(37) 



E ((cfJ-cfi)*^)^) 



k|<-Rn 







(a.s.) . 



Using Markov's inequality and (|35|) . we conclude for any bounded B and e > that 

(#,5^(0)) 



(38) 



which is summable by (|36p . Hence, by Borel-Cantelli, 



(39) 



E c<2>**, 







a.s. 



.zGyl-yl 
|z|>-R» 



Combining this with (|37p shows that the limit is the expected one (from the pointwise calcu- 
lation) also in this case, which yields the claim. □ 

Note that our argument is based on the potential selection of a subsequence of the original 
(deterministic) A- However, it also shows that the limit derived in Theorem [2] is the only 
point of accumulation along any deterministic subsequence of A. 

Remark 9. Randomisation of Meyer sets. A particularly relevant class of point sets 
in the theory of aperiodic order are Meyer sets, which are relatively dense sets A such that 
A — A is uniformly discrete. Such sets always have a diffraction measure with a non-trivial 
pure point part, with a relatively dense supporting set [SHE], despite the fact that Meyer sets 
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can have entrop}H. If modified by a family of random measures according to Theorem [2l the 
resulting diffraction still shows the original diffraction with its non-trivial pure point compo- 
nent, modulated by the function |Eq(1?)| 2 , in addition to the diffuse background originating 
from the added randomness. 

Let us look at consequences of Theorem [2] in terms of some examples. 

Example 6. Deterministic clusters. Let 5 C l rf be a finite point set, and consider 
Q = 5s = ^2 x& s^ x - Clearly, this completely deterministic case gives Eq(|J?|) = Eq(J?) = 5s 
and Kq({2*{2) = 5s * 5s, so that Theorem [2] gives 7^ = f<5g * 5s) * 7^ and 7^ = \5s\ 2 • 74, 
which is always true (rather than almost always) in this case. A particularly simple instance 
of this emerges from S = {a}, which effectively means a global translation by a. This leads 
to the relations 7^ = 7^ and 7^ = 7a> as it must. 

Example 7. Random weight model. Here, we consider J? = H5q, where H is a complex- 
valued random variable with a law \x that satisfies E^dffj 2 ) < 00 (hence also E M (|ff|) < 00). 
Clearly, this gives E Q (f2) = E M (iT) 5 and E Q {Q*Q) = E M (|#| 2 ) <5 , so that Theorem [2] results 
in the diffraction formula 

lf> = |E M (#)| 2 • j A + dens(.l) (E M (|tf| 2 ) - \E,(H)\ 2 ) • A (a.s.). 

The autocorrelation is clear from Theorem [2j 

Remark 10. Interpretation as particle gas. A widely used special case of Example [7J 
is the random occupation model, or '/l-gas'. Here, Q may take the value 5q (with probability 
p, for 'occupied') or (with probability 1 — p, for 'empty'). This gives the diffraction 

1a Q) = P 2 • lA + dens(yl) ■ p(l - p) • A (a.s.), 

which was derived in a similar setting in [J], and later generalised to Bernoulli and Markov 
systems [3], to systems with finite range Gibbs measures [11], and beyond [lOlHI]. 

The results of Examples [6] and [7] can be extended in many ways, some of which will be met 
later on. One further possibility consists in replacing a point by a 'profile', as described by an 
integrable function, or by a finite collection of such profiles, which could represent different 
types of atoms. The corresponding formulas for the autocorrelation and the diffraction are 
then easy analogues of the ones given so far. 

Example 8. Random displacement model. Consider the random measure Q = 5 X , 
where X is an Revalued random variable with law v. If A C M d is a Borel set, one has 

(E Q (f2))(A) = [ 5 x (A)du(x) = [ l A (x)du(x) = u(A), 

which shows that Eq(1?) = v. One also finds Eq(J?*J7) = v{W^)5q = 5q. Theorem [2] now 
results in the equations 

iA = * 7^ + dens(yl) (5q — v*u) (a.s.), 

ffi = \V\ 2 ■ % + dens(yl) (1 - \u\ 2 ) ■ A (a.s.), 



'The binary random tilings of Example [S] produce Meyer sets whenever b/a G Q. 
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which recovers Hof's result on the diffraction at high temperature [33]. 

In comparison, Hof's approach to the random displacement model [33] also uses the SLLN, 
but does not require the FLC property. Instead, he needs an ergodicity assumption on the 
underlying point set; compare also [32]. 

Remark 11. Extension of Theorem [2j The argument above is shown for FLC sets in 
a pointwise fashion, to make the result more transparent. However, it is clear that one does 
not need the FLC property itself. Indeed, it is sufficient to assume that the fixed point set 
A, relative to a chosen van Hove averaging sequence .4, possesses an autocorrelation that is 
a pure point measure of the form 7 = ^Z z& pf]{z)5 z with F a locally finite point set. An 
argument with local test functions will then still connect to the SLLN and thus avoid the 
need for ergodic assumptions on the underlying set A. 

With hindsight, it is rather clear that the formulas of Theorem [2] are robust, and should 
also hold for other point sets, such as those coming from a homogeneous Poisson process or 
from a model set based particle gas, as introduced in [7J. So, to complement our approach of 
this section, let us now consider ergodic point processes instead, meaning that also the set A 
becomes part of the random structure. 



5. Arbitrary dimensions: Point process approach 

Here, we are interested in the diffraction of certain random subsets of where we restrict 
ourselves to the situation that these subsets are self- averaging in a suitable way. This will 
be guaranteed by the ergodicity of the underlying stochastic process. One further benefit 
is that we are freed from details of the averaging sequence and the potential selection of 
subsequences thereof. It is convenient to start by putting ourselves in the context of random 
counting measures, which we now summarise in a way that is tailored to diffraction theory. 

5.1. Random measures and point processes. Let AA + denote the set of all locally finite, 
positive measures 4> on M d (where we mean to include the measure). That 4> is locally 
finite (some authors say <j> is boundedly finite or that <p is a Radon measure) means that, 
for all bounded Borel sets A, 4>(A) < 00. The space AA + is closed in the topology of vague 
convergence of measures (in fact, A4 + is a complete separable metric space, see |19t A 2.6]j). 
We let denote the u-algebra of Borel sets of AA + . The latter can be described as the 

cr-algebra of subsets of AA + generated by the requirement that, for all Borel sets A C M. d , the 
mapping <j) \— > <j)(A) is measurable; compare \37\ Chs. 1.1 and 1.2] for background. 

A random measure on M. d is a random variable from a probability space (0,J-,ir) into 
(A4 + , Let us write V{M + ) for the convex set of probability measures on AA + . The 

distribution of a random measure $ is the probability measure P = P$ £ V(AA + ), defined by 
P = 7r o <P~ l . In other words, P is the law of ^, written as C(&) = P. Note that, as soon as 
P is given or determined, one can usually ignore the underlying probability space. 



We refer to the second edition of this work throughout, which comes in two volumes 19, 20 . All results we 
need are also contained in the original one volume edition [18] . sometimes with a slightly different numbering. 
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For each t £ R , let T t denote the translation operator on IR , as defined by the mapping 
x i — ^ t + x. Clearly, one has TtT s = Ti +S , and the inverse of Tt is given by Tf 1 = T_ t . For 
functions / on the corresponding translation action is defined via Ttf = / ° T-t, so that 
T t f(x) = f(x — t). Similarly, for (ft £ let T x (ft ■= (ftoT- x be the image measure under the 
translation, so that (T x (f>)(A) = <ft(T- x (A)) = (ft(A — x) for any measurable subset A C R d , 



and (T x (j))(f) = J Rd f(y) d(T x (ft)(y) = J Rd f(x + z) d(ft(z) = (/>(T- X f) for functions. This means 



that there is a translation action of W 1 on A4 + . Finally, we also have a translation action on 
V{M + ), via (T X Q)(A) = Q(T^ X A) for any measurable A C M + . 

Our primary interest is in random counting measures. A measure (ft on M rf is called a 
counting measure if (ft(A) £ No for all bounded Borel sets A. These are positive, integer- 
valued measures of the form (ft = J2iei > where the index set / is (at most) countable and 
the support of is a locally finite subset of M d . The (positive) counting measures form a 
subset J\f + C .M + . We can repeat the above discussion of Ai + by restricting everything to 
J\f + . The vague topology on J\f + is the same as its topology inherited from Ai + , and its 
a-algebra of Borel sets XV+ consists of the intersections of the elements of E M + with N + . 
The concepts of the law of a random measure and the translation action by R d carry over. In 
particular, for x £ supp(0) with (f> £ AA + , T- X (f) corresponds to the counting measure obtained 
from (j) by translating its support so that x is shifted to the origin. 

A point process on M> d is a random variable <P from a probability space (0,^,^) into 
(Af + , Alternatively, a point process is a random measure for which 7r-almost all 9 6 <9 

are counting measures. Furthermore, it is called simple when, for 7r-almost all 8 £ 0, the 
atoms of (ft = &(0) have weight (or multiplicity) 1. 

In many instances, the point processes we are dealing with are simple. Whenever this 
happens, we feel free to identify point measures with their supports. In this case, the measures 
almost surely are Dirac combs of the form (ft = 5s with S C R d locally finite. Later on, we 
create compound processes in which an underlying point process is decorated with a random 
finite measure, and this will take us from M + to A4 + , which is also the reason why we 
introduced random measures above. 

For a random measure (or a point process) ^ with law P, the expectation measure ¥,p(<P) 
is defined by 



It is a measure on R which gives the expected mass (or number of points) that <P has in A. 
More precisely, in terms of the underlying probability space (0, J 7 , tt), one writes 



It is common in the probability literature (and we adopt this slight abuse of notation here, 
too) to suppresses the explicit dependence on (0,J-,ir) by simply writing <P for the general 
instance &(0) of the process <P. The latter is called stationary when its law P is translation 
invariant, which means that T t P = P o T_ t = P holds for all t £ M d . 



(40) 
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Remark 12. Intensity of a process. If P is stationary, we have T t Ep(#) = Ep(#) for 
all t G ~R d , whence Kp(<P) must be a multiple of Lebesgue measure (the latter being Haar 
measure on Consequently, 



where p G [0, oo] is usually called the intensity of P. Unfortunately, this term is already in 
use for the positive weights of Bragg peaks in diffraction theory. In the setting of simple 
point processes, p also has the meaning of a point density, averaged over all realisations of the 
process. In the ergodic case (see below for a definition), it is then almost surely the density 
of a given realisation in the usual sense. We thus often prefer to call p the point density of 
the simple point process or the density of the random measure. 

From now on, we always assume that p is finite. Let ^ : {0,J-,n) — > {X,E%) be a 
stationary random measure (where X = M + ) or point process {X = with law P. Then, 

(X , Ex , P) is a probability space with translation invariant probability measure P. In fact, 
we will usually simply assume that (X, Ex,P) is itself the probability space (or basic process) 
we are dealing with. In general, there will be several different spaces, and to keep track of 
the processes, we use the law of the basic process as an index. 

Let us recall that the random measure or point process & is called ergodic when (X, Ex,P) 
is ergodic as a dynamical system |62] under the translation action of R d , see below for more. 
In particular, we do not refer to strict ergodicity this way. 

5.2. Palm distribution and autocorrelation. Let P G V(M + ) be stationary with finite 
density < p < oo. The assumption p > is no restriction, since it is easy to see that a 
realisation of a stationary point process with intensity p = almost surely is the zero measure. 

Let lg, as usual, denote the characteristic function of the set B C M + , and choose a Borel 
set Act'' with < X(A) < oo. The Palm distribution P is the probability measure on N + 
that satisfies 



for any B G ZJtf+, compare [571 Ch. 4.4] or [381 Ch. 3] for background. Due to stationarity, 
Remark 1121 applies to ¥,p(<P(A)), whence the prefactor simplifies to (p\(A))~ 1 . Note that the 
sum under the integral runs over at most countably many points. Moreover, the definition 
does not depend on the actual choice of A. Intuitively, P describes the configuration <P as 
seen from a typical point in supp(^), with that point translated to the origin. Alternatively, 
in the case of simple point processes, one can think of P Q as the distribution of ^, conditioned 
on having a point measure at 0. This actually amounts to condition properly on an event of 
probability 0, which might need some further explanation. 

The first point of view can be made precise, at least in the ergodic case, as a limit, via 
sampling points in ^ over larger and larger balls, see |38[ Thm. 3.6.6] or [20, Prop. 13.4.1 and 
Prop. 13. 4. IV] as well as Eq. (|43p below. The second interpretation can be corroborated by 
conditioning # to have a point in a small ball around and then again taking a limit, see |20l 
Thm. 13. 3. IV]. In more precise terms, P would be called the Palm distribution with respect 



I P (#) = E P (#) = pX 



(41) 




xGAnsupp(^) 



30 



MICHAEL BAAKE, MATTHIAS BIRKNER, AND ROBERT V. MOODY 



to G M. d , compare |36|, Ch. 10] or |20} Ch. 13.1]. Since we will mostly be dealing with the 
stationary scenario, we refrain from spelling out the full name. 

There is an alternative approach to the Palm distribution, which also applies to the random 
measure case, compare [20], Chs. 13.1 and 13.2]. Let <P : 7r) — ► (M. + , Sm+) be a 

stationary random measure with law P and finite mean density p < oo. Then, the Palm 
distribution is the unique probability measure Pq on A4 + that satisfies 

(42) Ep( f g(x,<P)d<P(x)] = p [ [ g{x,T x i;)dP (i,)dx 

\JR d J JR d J M+ 

for all non-negative functions g on M d x A4 + for which J^ d j M+ g(x, <fi) d(j)(x) dP(4>) is finite. 
When dealing with point processes, all this reduces to M + by simply replacing M + with J\f + 
throughout Eq. (|42]), compare [20] Ch. 13.2 and Thm. 13.2.III]. 

If <P is an ergodic, stationary random measure, an application of the ergodic theorem implies 

(43) [ F(T^)d0(x) ^> pf F(V)dP (V) (a.s.), 

for any non-negative measurable function F : M + — > M, see [20, Prop. 13.4.1] or the proof of 
Thm. 3.6.6]. Here and below, we write B n for B n (0) and X(B n ) for vo\(B n (Q)). 



In the literature, the probability measure P Q is usually called the Palm distribution of P 
(with respect to 0), while the term Palm measure is also in use for the unnormalised version 
pP Q , a convention we adopt here. The first moment measure of the latter coincides with the 
autocorrelation measure of the underlying process and is denoted by j P . This is motivated 
by the following result on the autocorrelation 7p of a given realisation, which is somewhat 
implicit in the literature. Its importance in our present context was first emphasised by 
Gouere in [28]; see also [33] for complementary aspects. 

Theorem 3. Let <& be a stationary, ergodic, positive random measure with distribution P. 
Assume that P has finite density p, and that P has locally finite second moments in the sense 
that Kp(<P(A) 2 ) < oo for any bounded A C M d (this follows for instance from the condition 
Ep(^(_B r (x)) 2 ) < oo for all x £ IR d and some fixed radius r). Let <& n := ^\b u (o) denote the 

restriction of # to the centred ball of radius n. Then, the natural autocorrelation 7^ of 
defined via an averaging sequence of centred nested balls, almost surely exists and satisfies 

^ '= ^ vol(g w (0)) = iSc val(g w (0)) = pIp o = 7P ' 

where the limit refers to the vague topology on Ai + . Here, Ip is the first moment measure 
of the Palm distribution, 

I P (A) = [ &(A)dP (&), for A cR d Borel. 
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Proof. As test function, fix a non-negative continuous function g: M. d — > [0, oo) with compact 
support. With B c n := R d \B n , we have 

— ^— / g(x)d(<P n *¥ n )(x) = — ^— [ g(x - y) d<P(x) d$(y) 

A V-DnJ jRd A^nj JB n xB n 

= TiTT~\ I \ 9{x-y)6${y) - f g(x - y) d<P(y) J d$(x) 

= F g {T^ x <P)d<P(x) - R n {g) 

(note that both integrals inside the big brackets in the second line are finite because g has 
compact support), where (ft t— > F g {4>) = j Rd g(— z) dcft(z) defines a measurable function, and 
the remainder is given by 



Rn{g) = / / g{x - y) d${y) d$(x) 



Note that R n , which is a random measure, is precisely the difference between the elements of 
the two approximating sequences of random measures in the claim. In view of (|43p . it thus 
remains to show that linin^oo R n = almost surely. Choose k so that g(x) = for \x\ > k, 
and fix some e > 0. We then have, for n > k/e, 

R n (g) < jj^ / a»(BS n (x + B k )) d$(x) < f #(x + s*) d*(x) , 

where (j) ^ G((p) := <fi(Bk) is again measurable. Hence we obtain 

*\£>n) JB n *{r>n) A(±f(i_ e ) n J JB {1 _ E)n 

n ^°% (l n ~\d\ ||„|| „/ rifiTAA-DhTA — fl M „\<A 



(1 - (1 - e) d ) HfflU p / G{W) dP (#) = (1 - (1 - e) d ) Halloo pI P (B k ), 

almost surely by (|4"3"|) . Now take e\0 to conclude. □ 

Continuing with the hypotheses of Theorem O our assumptions guarantee that the second 
moment measure // 2 ) of P, defined on cylinder sets A x A' C M. d x M. d via fi^(A x A') = 
J^- + &(A) 4>(A') dP($>), is locally finite. This is a necessary and sufficient condition for the 
existence of the first moment measure of the Palm distribution (as a locally finite measure). 
In fact, in the stationary scenario, the autocorrelation of the process, denoted by "y p , satisfies 

(2) (2) 

7p = Mred' where fi re ^ is the so-called reduced second moment measure of P, and this, in turn, 
is the same as the intensity of the Palm measure. We offer a brief explanation of this (for 
more details, see [20] Prop. 13.2.VI] or [57J Ch. 4.5]). 

First, /^ed 1S obtained from by disintegration, via factoring out the translation invari- 
ance. More precisely, following [20], p red is the unique positive measure on Mr such that 

(44) / h(x,y)d^(x,y) = [ f h(u, u + v) d/j™(v) d\(u) , 

JR d xR d JR d JR d 
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for all (real) functions h £ C c (R d x R d ). When h = f <8> g is a product function (meaning that 
h(x,y) := f(x)g(y)), one finds the relation 

(45) P {2) (f®g) = i^(f*g) 

(2) 

via Fubini's theorem. Choosing g = f, it is clear that the measure A is positive definite. 
More generally, when dealing with complex-valued functions, one has to consider 

P (2 \I®g) = &f*g), 

which leads to some technical complications later on. Since we consider real- valued component 
processes only, we can stick to the simpler case of real-valued functions. 

The connection of the reduced second moment to the intensity measure of the Palm measure 
comes through applying (|42p to a function on M rf x A4 + defined by 

(46) (x,<j>) i ^ g(x) i T x h{y)dcp(y), 

JR d 

where g, h are arbitrary, but fixed, non-negative measurable functions onl <i . The left hand 
side of (H2]) then reads 



e p( / 9(x) I h(y - x) d${y) d${x) J = Ep\^J J g(x)h(y - x) d$(y) d$(x) 

g(x)h(y - x)d^ 2) (x,y) = \{g) ■ iifj A (h) , 
where we employed Fubini's theorem and (I44j) . while the right hand side reads 



P I j + 9(x) I (T x h)(y)d(T x 4>)(y)dPo(0)d\(x) 



g{x) / h(y) d<j>[y) dP (0) dA(x) = \{g) ■ P I Po (h) . 

M+ JR d 

Here, we used the notation of the intensity of the Palm measure for its first moment. Com- 
paring these two calculations gives 

(47) ^ d = p I Po=lp . 

(2) 

Remark 13. Equivalent definitions of fi y ie ' d . There are several different ways to define 
a reduced measure via disintegration. In particular, one could employ h(u, u ± v) as well as 
h(u ± v, u) in Eq. (|44|) . Using translation invariance of Lebesgue measure, this boils down to 
just two different possibilities, the one with h(u, u + v) introduced above and the one with 
h(u,u — v), which is used in [371 Prop. 1.60]. Observing the relation 

f*g = f*g 

together with /^( 2 ) = one can check that both versions define the same measure, as 

the process is restricted to positive (and thus real) random measures, so that no complex 
conjugation shows up in the 7-operation. Alternatively, one can use commutativity of the 
convolution together with the symmetry of p^ 2 \ which implies p^ 2 \f <S) g) = p^(g <8> /)• 



DIFFRACTION OF STOCHASTIC POINT SETS 



33 



Remark 14. Renewal process revisited. Consider a stationary renewal process on 
the real line, viewed as a random counting measure, with inter-arrival law g. The latter 
is assumed to be a probability measure on M + , with expectation J" R xdg(x) = 1. It is 
well known (compare [20\ Thm. 13.3.1 and Ex. 13.3(a)]) that the Palm distribution P , with 
respect to the origin, of (the law of) equals the law of 

(48) ^o+E%< 
where S = 0, 



S, 



\X 1 + ... + X i , if i > 1, 
[X i+1 + ... + X , ifi<-l, 

and (Xj) igZ are i.i.d. with law g. We see immediately from (j48]) that the first moment measure 
I p of P is given by formula (llip . see also Proposition [H and thus recover Theorem [T] as well 
as Remark [3] by specialising Theorem [3] to the case of a renewal process on the line. 

With the interpretation as a random counting measure, we are no longer restricted to laws 
g on M + . Indeed, when g is any probability measure on R with expectation 1 (which prevents 
the process from being recurrent), the random counting measure almost surely leads to the 
autocorrelation given in (jlip . and thus avoids the complications mentioned after Lemma [H 
see also [251 Ch. XI.9] for a systematic exposition, and [El Ex. 8.2(b)] for comparison. 

Remark 15. Bartlett spectrum. The diffraction measure 7 of a stationary random 
measure (with E(|<£(A)| 2 ) < 00 for all compact A C say) is closely related to the 
so-called Bartlett spectrum r := (c±S) of as follows; compare [HI Ch. 8.2]. 
Recall that the covariance measure is defined on cylinder sets via 

c {2) (AxA') = J 0(A)0(A')dP(0) - J &(A) dP($) J 0(A')dP(0) 

= ^ 2 \A x A') - p 2 (X 8A)(Ax A'), 

where p is the density of 0, compare [20, Eq. (9.5.12)]. Consequently, the relation between 

(2) (2) 
the reduced second moment measure p ied and the reduced covariance measure c re j of is 

(49) <ffl-/ffl-^- 

^ f<2~\ /n\ _____ 

Since 7 is the Fourier transform of fi red and r that of c^J,, Eq. (|39j) translates into 

(50) r = 7 -p 2 ,5 . 

From our perspective, the positive definite autocorrelation measure 7 is a slightly more natural 
and universal object than the inverse Fourier transform r of r, since 7 corresponds to a 
physically observable quantity, namely diffraction. Eq. (|50p gives r = 7 — p 2 X, which is no 
longer positive definite. Rather, it is tailored to situations where supports the only atom of 
7, as in the homogeneous Poisson process; compare Example [9] below and the discussion in 
[T9l p. 305], and see [191 Sec. 8.2] for further explicitly computable examples. 
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To formulate the standard Poisson process in this setting, let us start with an intuitive 
picture. Imagine independently putting single points on the sites of eL d C each with 
probability pe d , and imagine a process that arises from this construction in the limit e — ► 0. 
For a rigorous construction, one can start from a tiling of W 1 with translates of [0, l) d and 
then proceed, independently for each cell, as follows: Put a Poisson-(p) distributed number of 
points in the cell, with their locations independently and uniformly distributed over the given 
cell, see |57t Sec. 2.4.1] for details. Such a more elaborate approach is needed when d > 1, as 
there is no analogue of the renewal process we used for d = 1. 

Example 9. Homogeneous PoiSSON PROCESS. This process on R d , with (point) density p 
(compare Remark HJ), is a random counting measure <P (with distribution P) such that ^(^4) 
is Poisson-(/?A(^4))-distributed for any measurable A C M d and that the random variables 
&(Ai), . . . ,$(A m ) are independent for any collection of pairwise disjoint A\, . . . ,A m C M. d . 
With this setting, the expectation measure of the process is given by Ep(#) = pX. 

It is well-known that, under the Palm distribution, a Poisson process looks like the same 
Poisson process augmented by an additional point at 0, so that 

(51) P (B) = Jl B {$ + 6 )dP($), ioiBcN 

(alternatively, write C{<& + 5q) = P , or P * 5s = P ), by a theorem of Slivnyak, compare 
\57\ Ex. 4.3]. This is intuitively obvious from the approximation via independent coin flips 
on eZ d and the idea of obtaining the Palm distribution via conditioning on the presence of a 
point at 0. Here, this gives in Ip Q = 5o + Ip = 5 + pX. Since homogeneous Poisson processes 
are stationary and ergodic for the translation action of M d , we can now apply Theorem [3l 

Consequently, for almost all realisations <£> of a homogeneous Poisson process of density p, 
the autocorrelation measure and the diffraction measure are given by 

(52) 7 P = p5 + p 2 X and 7 P = p 2 5 + pX, 

by an application of Eq. Q . This also extends Example Q] to arbitrary finite values of the 
point density p; compare also [19j Ex. 8.2(a)]. 

Remark 16. Matern's hard-core process. One drawback of the Poisson process (of 
point density p > in with d > 2 say) for applications in physics is the missing uniform 
discreteness. The latter was introduced by Matern through a hard-core condition, realised via 
a local thinning process applied to each realisation, see [58] and references therein. Informally, 
for some fixed radius R > 0, his construction works as follows. Each point of a realisation of 
a Poisson process is equipped with an independent mark that is drawn uniformly at random 
from [0, 1] (technically, this is a marked Poisson process). Then, only those points x are kept 
for which there exists no point with a smaller mark in Bp{x). The autocorrelation of the 
modified process is still radially symmetric. Moreover, if R is the radius of the hard-core 
condition and B R = B R (0) as before, the autocorrelation for distances r > 2R is that of a 
homogeneous stationary Poisson process with a new, effective point density 

1 _ e -P vo K B R ) 

Peff = vol(B fl ) • 
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In fact, the new autocorrelation almost surely has the form 7 = p eff 5 + PeS^ ~ v ^ where v is 
a radially symmetric measure with support B 2R , as follows from |58^ Thm. 1]. In fact, v is 
absolutely continuous with density p^ s l B — g, where g is a smooth function on B 2R \ B R . 
This density has a jump (with sign change) at \x\ = R, but tends smoothly to as |x| /* 2R. 
The diffraction of (this version of) the Matern model is thus given by 

7 = PcfT^O + (Peff - h ) x ( a - s )> 

where h is a smooth (even analytic) function with h(k) = 0(|fc|~( d+1 )' 2 ) as \k\ — ► ±00. 
Indeed, one has 

/ R \ d/2 

lB n {k) = {]k\) J d/2(^\k\R), 

which is responsible for the above estimate via the exact asymptotic behaviour of the Bessel 
function J d j 2 for large arguments. The remaining contribution, using the explicit expression 
of |581 Eq. (3.2)] and the reduction of the radially symmetric Fourier transform to a one- 
dimensional Hankel transform, gives another term. It can also be computed explicitly, though 
the resulting formula is too lengthy to write it down here. Its decay is as for the previous 
term, by an application of the (refined) Paley- Wiener theorem. 

At this point, let us recall from [7] one process that is of particular interest in the study of 
aperiodic solids. Unlike the Poisson process and most of its siblings, it shows a substantial 
amount of point spectrum. It is related to Remark [TOj, but based on the cut and project 
formalism, for which we refer the reader to [44] . 

Example 10. Model set based particle gas. Let A C R d be a regular model set, for 
simplicity with internal space R m . Let C be the lattice in R d x R m that is needed for the cut 
and project scheme, with projection image L in R d . We denote the corresponding star map 
(from L into R m ) by *, so that (up to a translation) 

A = {x G L I x* G W}, 

where W C R m is the window of A. The latter is assumed to be a compact set with non-empty 
interior and a boundary of measure 0. This guarantees A to be a Meyer set. Let / now be a 
continuous function on W, and consider the weighted Dirac comb 

xEA 

which is pure point diffractive by the model set theorem [32} [551 \7\ with diffraction measure 

^ = E \K-k*)\ 2 h- 

fee M 

Here, M is the projection of the dual lattice C* into Mr, on which the star map is well-defined, 
too. It is known as the Fourier module of the model set A. 

Assume now that < f(y) < 1 on W, and define a family of independent binary random 
variables {U x ) x&A , each taking values or 1 with F({U X = 1}) = f(x*). The stochastic 
counterpart lo s of u> is 
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which can be interpreted as a particle gas on A. By Thm. 2 and Eq. (58)] , it almost surely 
has diffraction 

%j s = iZ + dens(yl) V X. 
Here, V = vol ( W ^ f w f(y) (l — f(y)) dy is the mean variance of the random variables, averaged 
over A, which is a consequence of the uniform distribution result for model sets [541 145] . One 
can also calculate the entropy of this system [7\. Moreover, it is not difficult to restrict the 
process to produce Meyer sets - one simply has to choose a function / that is 1 on a subset 
of W with non-empty interior. 

5.3. Compound processes. Let us now go one step further by adding random clusters to 
the picture. To this end, let a stationary, ergodic, point process be given, with law P, 
density p, and locally finite expectation measure Ep(#). This is called the centre process 
from now on. Moreover, let \P G -M-m be a positive random measure with law Q, subject 
to the condition that both its expected total mass, m := E Q (#(R d )) > 0, and the second 
moment of the total mass, E ( g((lP r (]R ci )) 2 ), are finite. This is the component process. We will 
also consider signed component processes & with values in M bd , in which case we assume 
that the second moment of the total variation measure is finite; see the appendix for some 
details on the required notions and modifications for signed measures. 

A combined cluster process, or cluster process for short, is a combination of a centre process 
and a component process of cluster type, and is obtained by replacing each point x G supp(#) 
by an independent copy of 1^, translated to that point x. We denote such a process by the 
pair (^p,1^q). As before, we restrict ourselves to finite clusters here. Formally, let !Zq,*p2, • • • 
be independent copies of \P (these are the individual clusters). When ^ = ^ 6 X ., we put 

i i 

and denote the resulting law by P c \ = R. Note that, when W = 5q is deterministic and 
concentrated to one point, we simply obtain £(#p,<?q) = £(&), and the cluster process 
coincides with the centre process. 

If ^ is a counting measure, the cluster process (^p,1^q) is again stationary and ergodic, 
and its expected density is given by mp, by [20, Prop. 12. 3. IX]. This property actually holds 
in larger generality, which we need later on. 

Proposition 3. Let be a stationary and ergodic point process with law P, finite density p 
and locally finite second moments. Let \P be a (signed) random measure with law Q, finite 
mean and finite second moment. Then, the combined cluster process, which is a random 
measure, is again ergodic. 

Sketch of proof . If the component process is a (positive) point process as well, this result 
is stated and proved in [20|, Sec. 12.3]. The necessary modifications for an extension to a 
(possibly signed) random measure as component process, which seem to be well-known but 
which we could not explicitly trace in the literature, are provided in the appendix. □ 

The second moment measures of the three processes are connected in a way that permits an 
explicit calculation of the autocorrelation 7^ in terms of 7 P and various expectation measures 
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of the component process with law Q. To employ this powerful connection, we recall another 
disintegration formula, this time for any random variable S of the cluster process, 

(53) Kr(S) = F,r(F,r(S | given the centres)) = Kp(Kq(S | given the centres)), 

which (with obvious meaning) follows from the standard theorems on conditional expectation. 
We are now in the position to use Eq. (|44p in conjunction with Theorem [3] and Eq. (|53|) 



(2) (2) 

to calculate fi^ = fx R , and thus the autocorrelation of almost all realisations of the cluster 
process, where we first concentrate on positive random measures. The extension to signed 
measures follows in Section [5.41 We first need some technical results. 

Lemma 8. Let A be Lebesgue measure on M. , as before, and \i a finite Borel measure. Then, 
one has \x * A = cA with c = id(M. d ). 

Proof. Let g be an arbitrary continuous function on M d , with compact support. For all x £ 
we have \(T_ x g) = (T x X)(g) = A(<?) due to translation invariance of A. The convolution (j,* A 
is well-defined as \x is finite while A is translation bounded [131 Prop. 1.13]. One thus has 

(fj,*\)(g) = / g(x + y) dA(y) d/x(x) = / \(T_ x g) dfi(x) 

JR d xR d JR d 

\(g)d»(x) = /i(M d )A( 5 ). 

Since g was arbitrary, the claim follows. □ 

Given a measure ^ G M + and a continuous function g on M. d with compact support 
(possibly complex- valued) , we define a new function g^ on ~K d via 

(54) g„(x) := (T x [i)(g), 

which is certainly measurable. It is easy to check that g^ satisfies 

(55) gjj, = fa 

Lemma 9. Let /i G Aiy and let 7 be a positive, translation bounded measure on M d . For 
arbitrary (possibly complex-valued) f,g € C C (R ), one has the identity 

(/i*/2* 7) (/*<?) = 7(/m*^u)- 
This identity also holds when both /i and 7 are signed measures. 

Proof. Let / and g be //-measurable functions such that / * g is a continuous function with 
compact support, which includes the situation of the claim. One then finds, with Fubini, that 

[»*fi*l)(f*9) = J J(J f(x + z + £)Mx))( /$(v-0d£(v))dA(0d7(*) 

{T zHf i)(f) (T_^)(s)dA(£)d 7 (*) 



+ 0&(-0dA(0d 7 (*) = 7(/m*5~m) 
where all integrals are over M d and (|55l) was used in the last step. □ 
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Specialising Lemma [9] to 7 = 5q gives the relation 

(56) {ii*fi)(f*g) = (/m*&)(0) = KfM, 

which simplifies our further discussion. 

Remark 17. Test functions for measures. Recall that two measures /i, v G .A/f are 
equal when fi(h) = v{h) for all h G C C (M ). When the measures are positive or signed (but 
not complex), real- valued functions suffice. In the latter case, it will be particularly helpful 
to restrict to functions of the form h = / * g with f,g£ C c (K d ). Since the space C c (IR d ) 
contains an approximate unit for convolution, the linear combinations of such functions are 
dense in C c (M. d ), so that they suffice to assess equality of measures. 

Lemma 10. Under our general assumptions on the component process, one has 

(E Q (#**))(f*g) = A(E Q (/^)) and 

(E Q (#)*E^))(/*$) = A(/ Eg( , )S ^j), 

where f,g € C C (R ), possibly complex-valued. 

Proof. Let / and g be chosen as in the previous proof, with complex-valued functions per- 
mitted. For the first claim, observing that each realisation of 'F is a finite measure, a direct 
calculation with Eq. (|56p gives 

(Eq(#* #))(/*$) = Eq((#* £)(/*#) = Eq(A(/„^)) = A(Eq(/^)), 

while the second identity simply is Eq. (|56p with ^ = Eq(^), which is a finite measure by 
assumption. □ 

Recall that the covariance of two real- valued random variables X and Y related to the law 
Q (using our general notation as explained above) is defined as 

(57) cov Q (X,Y) := E Q (X Y) - E Q (X) E Q (Y) , 
where the index Q highlights the reference to the underlying law Q. 

Proposition 4. Let (S, R) be a combined cluster process with stationary centre point pro- 
cess (&,P) and real component process (<P r , C}), both with the usual assumptions on means 
and second moments as used above. Then, S is locally square integrable in the sense that 
E^[(^(i?)) 2 ] < 00 for any bounded Borel set B, and we have the reduction formula 

V ( R(f®g) = (/e q (*) ®5 Eq (^)) + P X { cov q(U'9#)), 

where p is the density of the centre process, f, g are continuous with compact support, and the 
covariance is defined as in (|57p . 

Proof. In order to check that E^j [(S"(i?)) 2 ] < 00 for bounded, Borel measurable B C one 
can trace through the steps below, replacing / and g by 1 B (the corresponding integrals then 
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involve only positive terms and are finite by Fubini's theorem). In general, by assumption 
and the disintegration formula (j53[) . one finds 

^\f®9) = [ B(f)S(g)dR(S) 
JM+ 

= [ E Q ( Yl ^(r- a; /)^(r_ 2/5 ))dP(#), 

x,j/Gsupp(s>) 

where \P X denotes the random measure at centre x. Since \P X and Hs y are independent for 
x ^ y, the double sum over the support is split into a sum over the diagonal (x = y) and a 
sum over all remaining terms {x 7^ y). Using the linearity of the expectation operator, the 
integrand can now be rewritten as a sum over two contributions, namely 

E Q {V(T- X f)) Eq (#(T-v9)) and 
(Eg [HT-xf) HT- X g)) - Eg [9(T-xf)) Eq (HT-xg)) ) • 

X 

Inserting the first term into the previous calculation leads to the contribution 

^p ] { e q(U) ®Eq(s^)) = (/eq(*) ®^E (*))> 
while the second results in 

E P (#)(cov Q (/^,^)) = pA(cov Q (/^, 5 ^)), 
where the last step follows from the stationarity of (<P,P). □ 

Theorem 4. Lei $ be a stationary and ergodic point process with law P , finite density p 
and locally finite second moments. Let \P be a random measure with law Q, finite expectation 
measure and finite second moments. If (S, R) denotes the combined cluster process built from 
the centre process (&,P) and the component process (^,(5), it is also stationary and ergodic. 
Moreover, the autocorrelation of the combined process satisfies 

lR = (E Q ($)*E^(¥)) * 7p + p(Eq(^)-Eq(^)*E^)), 
and this is almost surely the natural autocorrelation of a given realisation of the cluster process. 

Proof. Choose two measurable functions / and g such that f*g£ C c (M. d ). Then, in line with 
Remark [T7] and Eq. (|45p , one finds via Proposition [J] that 

7p(/*<?) =VR > (f®g) = (/e (*) ® 9e q (*)) + P x { cov q(.Ui9*)) 

= 7p(/ Eq (*) * %^ } ) + P (Eq(^* - E Q (V) * E^)) (f*g), 

where Ep(<P) = pX due to stationarity of (#,P). The second step makes use of Lemma [TUl 
The formula for the autocorrelation now follows from the observation that 

7p(/e q (*)*^5)) = (Eq(^)*Eq(^)* 7p )(/*5), 

which is an application of Lemma [9j The remaining claims are clear due to the assumed 
ergodicity, via an application of Proposition [3l □ 
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An application of the convolution theorem gives the following consequence, where also the 
identity Eq(<P') = Eq(^) was used to highlight the structure of the result. 

Corollary 1. Under the assumption of Theorem^ the diffraction measure of the combined 
cluster process is given by 

% = |e q (^)| 2 -7p + p(iEo(^I 2 )-|iEqWI 2 )a 

which is then almost surely also the diffraction measure of a given realisation. □ 

The result parallels our previous formulas, as was to be expected. Nevertheless, it does 
not follow from Theorem [2] in general, because realisations of stationary point processes in 
M rf generically fail to be FLC sets. Before we discuss possible generalisations beyond the case 
of positive random measures, let us look at some examples. 

Example 11. PoiSSON cluster process. An important special case emerges when the 
centre process is the homogeneous Poisson process of Example El with point density p. Let 
7p and 7p be the corresponding measures. If we couple a cluster component process \P to 
it, with law Q and m := Eg(lP r )(]R ci ) its expected number of points, our general formula for 
the compound process (<Pp,<Pq) applies. With Lemma the convolution formula can be 
simplified, and the result reads as follows. 

For almost all realisations of a Poisson cluster process (&P, &q), the natural autocorrelation 
measure exists and is given by 

1r = 1p,q = ijnpfX + P E Q (M) , 

where Kq^*^) is a finite positive measure (of expected total mass > m 2 ), due to our general 
assumption that Eq^^I^)) 2 ) is finite. Consequently, the diffraction measure is almost 
surely given by 

7 K = (mp) 2 5 + p(E Q {M)) -A, 

where (Eq is a uniformly continuous Radon-Nikodym density for Lebesgue measure. 
These formulas include the case of deterministic clusters; compare Example [6l 

Remark 18. Random displacement of Poisson processes. An interesting pair of 
processes is the combination of the homogeneous Poisson process from Example [9] with Hof's 
random displacement model from Example O A simple calculation shows that 

Ir = 7p = 7p and Tp = 7 P 

in this case (and, in fact, R and P have the same law here). From a physical point of view, 
this is in line with the behaviour of an ideal gas at high temperatures. When the Poisson 
process is a good model for the gas, and random displacement one for the disorder due to 
high temperature, compare the discussion in [33J , the combination should still be an ideal gas 
- and this is precisely what happens, as reflected by the two identities. 

Remark 19. Particle gas cluster process. It is clear that the particle gas of Ex- 
ample [TU] satisfies all requirements for a centre process, so that we can apply the cluster 
process machinery to it, too. This produces physically interesting and relevant examples with 
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a substantial amount of point spectrum. This observation remains true for more compli- 
cated particle gas models with interactions, under certain conditions on the potential of the 
underlying Gibbs measure, say; compare [11] for further details and examples. 

Example 12. Neyman-Scott processes. Let K be a non-negative random integer with 
law C{K) = fj,, mean m := E^K) and finite second moment, E^(K 2 ) < oo. Now, let 
Y%,Y2, ... be a family of M. -valued i.i.d. random variables with common distribution u, and 
independent of K. Define the cluster distribution via 



K 

v.- 



9 := J> 

3=1 

i.e., a cluster has a random size K, while the positions of its atoms are independently drawn 
from the probability distribution v. The induced distribution for 9 is again called Q. With 
a calculation similar to the one in Example El one finds 
K K r 

Eq(¥)(A) = E Q (j2U{Xi^ = E M (£ / lA(xi)du(xii) = E„(K-v{A)) = mv(A) 
i=l ' i=l J Rd 

for A C M d Borel, so that Eq(W) = rnu and Eq(\P) * Eq(iF) = m 2 {v * v). Moreover, one has 

K 

Eq(V*V)(A) = E q (^ l A (X k -X e )) = mS (A)+M ll [K(K-l))(u*v)(A) 1 

k,e=i 

which gives Eq(\P* 9) = m5o + E^fK^K — 1))(^ * u), so that the general formulas from 
Theorem [2] can now be applied again. Note that E^KiK — 1)) = E^K 2 ) — m. 

If the centre process is once more the homogeneous Poisson process with mean (point) 
density p, Lemma [8] gives similar simplifications as in Example [TT1 Consequently, for the 
resulting law R, the autocorrelation is almost surely given by 

7 fi = (rap) 2 A + mpS + p(E^(K 2 ) - m) (y * v) , 

whence the corresponding diffraction measure is given by 

7 H = (mp) 2 5 + p(m + (E^K 2 )-m)\u\ 2 ) A, 

which is an interesting extension of the Poisson process; compare [191 Ex. 8.2(f)] for a circu- 
larly symmetric case in M 2 . 

5.4. Autocorrelation for signed (ergodic) processes. It is intuitively clear that the 
results of this section are not really restricted to point processes or positive measures for 
the clusters. Here, we sketch how they can be adapted to the situation of signed random 
measures. Consider a stationary, possibly signed, random measure H/ (with law Q and 'finite 
second moments', meaning that Eg((|lP r |(A)) 2 ) < oo holds for any bounded A C R rf ), with 
second moment measure p^ 2 \ defined as before via bounded / of compact support as 



/ 



R d xR d 



f(x, y) dfi^ ( x , y) = Eq( I f(x, y) d*(x) d*(yj) . 
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The reduced second moment measure n rcd on M d with the property 

(58) &U*9) = » {2) U®9) 

is defined in complete analogy to the positive case. The analogue of Theorem [3] is: 

Theorem 5. Let <P be a stationary and ergodic, random, signed measure with distribution 
P. Assume that <P has finite second moments in the sense that Ep((|#| (j4)) 2 ) < oo for any 
bounded measurable set A C M d (which follows, for example, from Kp^\^>\(B r (x))Y) < oo 
for all x S M rf and some open ball B r ). Let & n := ^|_b„ denote the restriction of <& to the 
ball of radius n around 0. Then, the natural autocorrelation of <P, which is defined with an 
averaging sequence of nested, centred balls, almost surely exists and satisfies 

'P := llm \/R \ = llm T77TT = Mred = 7P> 

where the limit refers to the vague topology on J\f . Here, fj%~\ is the reduced second moment 
measure of P according to (|58p . 

Proof. The proof is a variation of that of Theorem O Fix a continuous function h : M rf — > M. 
with compact support. We have to check that 

(59) A(fc) ( ^*^ )(/l) ^ {k) (a - S - } - 

Let ^ be an ergodic, random, signed measure as above and F an ergodic random function 
on the latter with the property that 

(60) E P ( / |F(x)|d|<2>|(x)) < oo 

J A 

for any bounded measurable A C M. d . We can then define an additive covariant spatial process 
Xa in the sense of |46j . indexed by bounded measurable subsets A, via 



X A := / F(x)6$(x). 
J A 



Note that ergodicity of ^ and F implies that (XA,^ d ) is again ergodic, meaning that the 
shift-invariant er-field is trivial. Now, Cor. 4.9] yields 

Applying this to <P as in the theorem, together with F(x) := L d h(x — y) d<P(y), yields 

If 1 ~ 

lim T77fT/ ^)d<£(x) = lim — — (£ n *<Z>)(/i) 

n-+ooA(±/ n ) J Bn n^oo X(B n ) 

= E ^(t777t/ / Kx-y)d$(y)d${x) 



1 



A (-Bl) JR d xR d 



1b 1 (x) h(x - y)dfi i2) (x,y) 



1 



A(5i 



l Bl (x)/»(«)dM2(«)da; = / hdv™ 
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almost surely, which is almost the claim. The difference between <P n * $ and <fr n * <P n can be 
treated as in the proof of Theorem [31 □ 

Combining Proposition [5] and Theorem [5j and observing that the calculations in the proof 
of Proposition 0] carry over literally to the signed case, we obtain 

Corollary 2. The statements of Theorem^ and Corollary{l\ remain true for cluster processes 
with signed clusters. □ 

Example 13. Signed PoiSSON process. If we combine the homogeneous Poisson process 
of Example [9] with the random weight model of Example and chose weights 1 and — 1 with 
equal probability, Corollary [2] implies the almost sure diffraction 

7 = pA. 

In particular, one has 7 = A for density p = 1, which makes this signed Poisson point set, 
on the level of the 2-point correlations, indistinguishable from the signed Bernoulli sequence 
(or process) on Z rf . This is remarkable in view of the rather different geometric structure and 
demonstrates the intrinsic difficulty of the corresponding inverse problem. 

5.5. Equilibria of critical branching Brownian motions in d > 3. Consider a system 
of particles performing independent Brownian motions in M. d , d > 3 (for ease of comparison 
with the cited literature, we assume that the variance parameter is a 2 = 2). 

Additionally, each particle, after an exponentially distributed lifetime with parameter V, 
either doubles or dies, where each possibility occurs with probability 1/2. In the situation of 
a birth event, the daughter particles appear at the position of the mother. Note that if we 
start with a finite number of particles, the expected number of particles is preserved for all 
time, as the expected number of offspring equals 1. This is what 'critical' in the name refers 
to. Imagine we start such a system from a homogeneous Poisson process of density p, denote 
by <&t the random configuration observed at time t > 0, and its distribution by Pt- Here, Pt 
is stationary with density p, see |31] and the references given there for background. 

It follows from [31[ Thm. 2.3] that the first moment measure of the Palm distribution of 
Pt is given by 

(61) I(P t)o = S + (p + f t )X, 

where ^ 

ft(x) = V / / p s {0,y)p s (y,x)dyds = — / p u (0,x)du, 

JO JR d 1 Jo 

with pt(x, y) = (A-Kt)~ d l 2 exp ( — \x — y\ 2 / (4t)) the d-dimensional Brownian transition density 
(with variance parameter 2). As explained in [31] , there is a genealogical interpretation behind 
f)61|) : In view of the interpretation of the Palm distribution as the configuration around a 
typical individual, 5q is the contribution of this individual, ft A that from its relatives in the 
family decomposition of the branching process, and p\ is the contribution from unrelated 
individuals. 

Furthermore, by |3H Thm. 2.2], Pt converges (vaguely) towards P^, which is the unique, 
ergodic, equilibrium distribution of density p (cf. [15] for uniqueness), and the limit t — ► 00 
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can be taken in (|6ip to obtain 

J (Pc=c)0 = 6 + (P+/oo)A. 



where 



v r°° v r(^Z) 1 

foo(x) = ~ / Pu (0,x)du = - [ ' ' 



2 J ^">^~" 2 4^ 2 |x| d - 2 
is (up to the prefactor V/2) the Green function of Brownian motion. Thus, using Lemma [2j 
we have 

Theorem 6. Let (fr^ be a realisation of the critical branching Brownian motion, from the 
equilibrium distribution P^. The autocorrelation is then almost surely given by 

7 = p5 + p(p + /oo)A, 

while 

is the corresponding diffraction measure. □ 

Remark 20. Extension of Theorem [6] One can also consider the scenario where, instead 
of Brownian motion, particles move during their lifetime according to a symmetric, stable 
process of index a 6 (0, 2] in R d (a = 2 corresponds to Brownian motion). Such processes 
have discontinuous paths, and their transition density pf\x, ,y — x) satisfies 



I e ihx p[ a) (0,x)dx = exp(-t\k\ a ) 



(in general, no explicit form of pf^ is known). By [31, Thm. 2.2], non-trivial equilibria exist 
if the spatial dimension d satisfies d > a. In this case, a reasoning analogous to that above 
yields the following: The autocorrelation of a realisation of the equilibrium of a system 
of critical, branching, symmetric a-stable processes (with density p) is almost surely given by 

7 = p 5 + p(p + fW)\, 

where 

^ V f 00 (n)/ , , V T((d-a)/2) 1 
•'oo v ) 2 y o F U v , > 2 2«7rd/2r(a/2) \x\ d ~ a 

(for the form of the Green function of the symmetric a-stable process, see |14l Ex. 1.7]). 
Hence, the diffraction measure is almost surely given by 

7 = P 2 <5n + pfl + — -, — ; — 7TT- ^A, 
' H P V 2 (2Tr) a \k\ a J 

by another application of Lemma [2] 

Note that, due to the independence properties of the branching mechanism, these equilibria 
can also be considered as Poisson cluster processes. In contrast to the scenario considered 
above, clusters in are infinite, and the spatial correlation decays only algebraically (with- 
out being integrable). 
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6. Outlook 



This article demonstrates that various aspects of mathematical diffraction theory for ran- 
dom point sets and measures can be approached systematically with methods from point 
process theory, as was originally suggested in [28]. At the same time, the approach is suf- 
ficiently concrete to allow for many explicitly computable examples, several of which were 
presented above. They comprise many formulas from the somewhat scattered literature on 
this subject in a unified setting. There are, of course, many more examples, but we hope that 
the probabilistic platform advertised here will prove useful for them as well. 

The next step in this development needs to consider point processes and random measures 
with interactions, such as those governed by Gibbs measures. First steps are contained in 
[331 HI [281 US HD (ini (ZD EE] and indicate that both qualitative and quantitative results are 
possible, though some further development of the theory is needed. 

A continuation along this path would also make the results more suitable for real applica- 
tions in physics and crystallography, though it is largely unclear at the moment what surprises 
the corresponding inverse problem might have to offer here. 

Appendix: Ergodicity for cluster processes with signed random measures 

Let A4 = A4(M rf ) be the space of (locally finite) real or signed measures on M , equipped 
with the topology of vague convergence, with M. + = A4 + (K d ) denoting the subspace of 
positive measures. Let Um denote the Borel cx-algebra of M d . Note that the latter is also 
generated by the mappings Ai 9 i— ► fi(A), for bounded and measurable sets A C Recall 
that any \i £ A4 admits a unique Hahn-Jordan decomposition 



The mappings p, \— > n+ and fj, 1— > jj,- are IT^-measurable. We write \n\ := /i+ + jX- £ M + 
for the total variation measure of fi. A random signed measure <P is a random variable with 
values in (A4, Em)- I n the context of signed random measures, it is convenient to work with 
the characteristic functional 



which is defined for any h: M d — ► K that is bounded and measurable with compact support. 
Here and below, we suppress M rf as the integration region. In analogy to the Laplace functional 
for positive random measures, the distribution of is determined by 92^. 

Here, we are interested in signed cluster processes: Let tf> be a stationary counting process 
with finite intensity p, and \Pj (with j £ N) independent (and independent from <£), identically 
distributed, random, signed measures such that E[|!Pi|] is a finite measure. Then, given a 
realisation <P = Y^j $Xj, where Xj are the positions of the atoms of <P (in some enumeration), 
the cluster process is defined as 



with p + , /i_ £ M. + mutually singular. 



(62) 



tp$(h) := E exp(a'f hd<Z>) 



(63) 
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Note that for any bounded B C R rf , 
E[|S(B)|] < E \J2 -Xj)] = p [ [ dE[\^\]dx = p(E[\^\]*\)(B) < oo , 

L "~ J JR d JB-x 

so that (1631) is indeed well-defined. 



Lemma 11. Let \P be a signed random measure on ~R d . The following are equivalent: 

(1) & is ergodic. 

(2) For any U,V G 27m» 

lim T7TTT / T X F) - P(# G U) P(lf G V)) dx = . 

(3) For any o, h: R d — > R measurable with compact support, 

1 



lim 

n^oo X(B n ) 



{<pv(g + T x h) - <p#(g) <p#{h)) dx = . 



Furthermore, it suffices to restrict to U,V to a semiring which generates Uj^ in (2), and it 
suffices to restrict to continuous g,h with compact support in (3). 

Proof. This is a straightforward adaptation of the proofs of Propositions 12. 3. Ill and 12. 3. VI 
and Lemma 12. 3. II of [20] to the signed case. □ 



The following result is an analogue \20\ Prop. 12. 3. IX] for the signed measure case. Since 
we have not been able to find a proof in the literature, we provide a sketch. 

Proposition 5. Let<P, and E := ^2jT x .^Fj be as above. If<P is ergodic, then S is ergodic 
as well. 

Sketch of proof. We verify condition (3) from Lemma fTTT Observe that for any / : R d — ► R 
with compact support and any e > 0, we can find R < oo such that 



(64) 



.j:\Xj\>R 



>£]<£. 



To check let R' be large enough so that supp(/) C [—R r , R'] d , and note that for R > R', 
the left-hand side of ([Ml) is bounded by 



< 



■E 



» j ■ l-Jvj |oo: 

The expectation on the right-hand side above equals 



£ \T x W 3 \([-R',R'] d ) 

lj:\X j \ cc >R 



P 



y)dE[\V 1 \](y)dx 



l d \[-R,R] d R d 



< p(2R') d E[\^ 1 \] (R d \ [-(R - R'), (R - R')] d ), 
which converges to as R — > oo because E is a finite measure. 



Let g, h : 
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Decompose 
j\g + T x h) d a 



E 



exp(ij gdS] 



H{$) := E exp(ifhdS) 



(p 



E / T ^9^ + E / 

+ E / Txj+xhdVj + 



j:X^[-i?,fl] d -a: 



T X . +x hd^, 



and choose i? so large that (|M|) is fulfilled for f = g and f = h. Recall that, for any real- valued 
random variables X, Y with P(|Y| > e) < e, we have 

E e i ( x+y ) - E e lX 



< E 


e i(X+Y) _ e iX 


< E 


'|e^||e iY -l|" 











< E- 



Y >e)< 2e. 



For A C write := Ylj- x^A^x^o f° r ^ ne ran dom measure which consists of clusters 
with centres in A. For x £ M. d \ [-2R, 2R] d , we then have 

E [exp {ij(g + T x h) d -)] - E [G(&)H(T a! &)] 

< 



E [exp (i f (g + T x h) dE)] -E exp (if g dS { _ RR]d + if T x hdS { _ RR]d _ x 



+ 



E 



E 



exp (if gdS { _ RJi]d + if T x hdS { _ RR]d ^ x 



P 



E [G{<P)H(T X <P)\ 



The first term on the right-hand side is bounded by 2e. Observing that the conditional 
expectation in the second term is in fact a product because clusters with centres in disjoint 
regions are (conditionally) independent, we can bound the second term from above by 



E 



E 
+ 

< E 



exp (if gdS { _ RR]d 



P 



E 



exp [if T x hdS { _ RR]d _ x 



P 



H(T X <P) 



E 



E 



exp (if gdE { _ RR]d 



P 



G(<P)\ H(T X P>) 



exp (ifT x hdS[_ RRi d_ x ) - exp (i f T x hdS 



+ E 



exp (ifgd»! R 



R}< 



exp (if g dE) 



which is not more than 2e. 

Thus, using the relation E [E [exp (i f(g + T x h) dS") I #]] = tp s (g + T x h) together with 
E [G(&)] = ip B (g) and E [H{P>)\ = E [H(T x <P)j = ip s (h), we obtain " 

1 



lim sup 



(65) 



A(-B n 
< lim sup 



(<p#(g + T x h) - ip^(g) ifv(h)) dx 
1 



X(B r , 



E [G(<P)H(T X <P)] - E [G(<P)} E [H($)] ) dx 



4e = 4e 



by ergodicity of P> (in order to deduce this literally from statement (2) in Lemma QTJ one can 
for instance discretise the support of g and h and approximate G(<P), H(<P) with functions 
depending only on the random vector (^(cj))i<j<jv, where {q | 1 < i < N} is a collection of 
disjoint (small) cubes). Finally, take e — > to conclude. □ 
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